A&A  505,  385-404  (2009) 

DOI:  10.1051/0004-6361/200912041 

©  ESO  2009 


Astronomy 

Astrophysics 


The  large  quasar  reference  frame  (LQRF) 

An  optical  representation  of  the  ICRS 

A.  H.  Andrei^’^’^,  J.  Souchay^,  N.  Zacharias"^,  R.  L.  Smart^,  R.  Vieira  Martins'’^,  D.  N.  da  Silva  Neto^’^, 
J.  I.  B.  Camargo^,  M.  Assafin^,  C.  Earache^,  S.  Bouquillon^,  J.  L.  Penna^,  and  F.  Taris^ 

'  Observatorio  Nacional/MCT,  R.  Gal.  Jose  Cristino  77,  CEP2092 1-400,  RJ,  Brasil 
e-mail:  oatl@on.br 

^  Observatorio  do  Valongo/UFRJ,  Ladeira  Pedro  Antonio  43,  CEP20080-090,  RJ,  Brasil 
^  Observatoire  de  Paris/SYRTE,  61  Avenue  de  I’Observatoire,  75014  Paris,  France 
US  Naval  Observatory,  3450  Massachusetts  Av.  NW,  Washington,  DC  20392,  USA 
^  INAF/Osservatorio  Astronomico  di  Torino,  Strada  Osservatorio  20,  10025  Pino  Torinese,  Italy 
®  Centro  Universitario  Estadual  da  Zona  Oeste,  Av.  Manuel  Caldeira  de  Alvarenga  1203,  CEP23070-200,  RJ,  Brasil 

Received  12  March  2009  /  Accepted  20  May  2009 


ABSTRACT 

Context.  The  large  number  and  all-sky  distribution  of  quasars  from  different  surveys,  along  with  their  presence  in  large,  deep  astro- 
metric  catalogs,  enables  us  to  build  of  an  optical  materialization  of  the  International  Celestial  Reference  System  (ICRS)  following  its 
defining  principles.  Namely:  that  it  is  kinematically  non-rotating  with  respect  to  the  ensemble  of  distant  extragalactic  objects;  aligned 
with  the  mean  equator  and  dynamical  equinox  of  J2000;  and  realized  hy  a  list  of  adopted  coordinates  of  extragalatic  sources. 

Aims.  The  Large  Quasar  Reference  Frame  (LQRF)  was  built  with  the  care  of  avoiding  incorrect  matches  of  its  constituents  quasars, 
homogenizing  the  astrometry  from  the  different  catalogs  and  lists  in  which  the  constituent  quasars  are  gathered,  and  attaining  a 
milli-arcsec  global  alignment  with  the  International  Celestial  Reference  Frame  (ICRF),  as  well  as  typical  individual  source  position 
accuracies  higher  than  100  milli-arcsec*. 

Methods.  Starting  from  the  updated  and  presumably  complete  Large  Quasar  Astrometric  Catalog  (LQAC)  list  of  QSOs,  the  initial 
optical  positions  of  those  quasars  are  found  in  the  USNO  Bl.O  and  GSC2.3  catalogs,  and  from  the  SDSS  Data  Release  5.  The  initial 
positions  are  next  placed  onto  UCAC2-hased  reference  frames,  This  is  followed  by  an  alignment  with  the  ICRF,  to  which  were  added 
the  most  precise  sources  from  the  VLBA  calibrator  list  and  the  VLA  calibrator  list  -  when  reliable  optical  counterparts  exist.  Finally, 
the  LQRF  axes  are  inspected  through  spherical  harmonics,  to  define  right  ascension,  declination  and  magnitude  terms. 

Results.  The  LQRF  contains  100,165  quasars,  well  represented  accross  the  sky,  from  -83.5  to  -1-88.5°  in  declination,  being  10  arcmin 
the  average  distance  between  adjacent  elements.  The  global  alignment  with  the  ICRF  is  1.5  mas,  and  the  individual  position  accura¬ 
cies  are  represented  by  a  Poisson  distribution  that  peaks  at  139  mas  in  right  ascension  and  130  mas  in  declination.  As  a  by-product, 
significant  equatorial  corrections  are  found  for  all  the  catalogs  used  (apart  from  the  SDSS  DR5),  an  empirical  magnitude  correction 
can  be  discussed  for  the  GSC2.3  intermediate  and  faint  regimes,  both  the  2MASS  and  the  preliminary  northernmost  UCAC2  positions 
are  shown  of  astrometry  consistent  with  the  UCAC2  main  catalog,  and  the  harmonic  terms  are  found  to  be  always  small. 
Conclusions.  The  LQRF  contains  J2000  referred  equatorial  coordinates,  and  is  complemented  by  redshift  and  photometry  informa¬ 
tion  from  the  LQAC.  It  is  designed  to  be  an  astrometric  frame,  but  it  is  also  the  basis  for  the  GAIA  mission  initial  quasars’  list,  and 
can  be  used  as  a  test  bench  for  quasars’  space  distribution  and  luminosity  function  studies.  The  LQRF  is  meant  to  be  updated  when 
new  quasar  identifications  and  newer  versions  of  the  astrometric  frames  used  are  realized.  In  the  later  case,  it  can  itself  be  used  to 
examine  the  relations  between  those  frames. 

Key  words,  catalogs  -  reference  systems  -  quasars:  general  -  methods:  data  analysis  -  astrometry 


1.  Introduction 

The  establishment  of  the  ICRS  as  the  celestial  reference  sys¬ 
tem  (Arias  et  al.  1995;  Feissel  &  Mignard  1998)  answers  not 
only  the  desiderata  in  terms  of  the  choice  of  the  physical  model 
and  coherence,  but  also  of  practicality  when  confronted  with 
the  earlier  dynamical  model.  This  is  especially  true  since  the 
system  is  materialized  by  the  direction  of  quasi-inertial  grid 
points,  which  are  given  by  distant,  extragalactic  sources,  in¬ 
stead  of  nearby,  moving  stars.  However,  for  a  long  time  after 
the  establishment  of  the  ICRS,  the  faintness  of  the  extragalactic 


*  The  LQRF  is  only  available  in  electronic  form  at  the  CDS  via 
anonymous  ftp  to  cdsarc.u-strasbg.fr  (130.79.128.5)  or  via 
http://cdsweb.u-strasbg .  fr/cgi-bin/qcat?J/A-i-A/505/385 


sources  prevented  a  precise  determination  of  optical  positions 
for  a  large  number  of  them.  The  VLBI  radio  technique  answered 
this  need,  by  providing  precise  positions.  However,  such  precise 
positions  could  only  be  derived  for  a  limited  number  of  strong 
and  stable  sources,  and  the  ones  with  richer  observational  history 
have  come  to  define  the  ICRF  (Ma  et  al.  1998).  Yet  even  this 
small  number  of  sources  enabled  the  VLBI  to  make  important 
contributions  in  defining  both  geophysics  and  radio  astronomy 
standards. 

In  contrast,  for  extragalactic  sources  on  the  optical  domain, 
poor  astrometry  and  a  small  number  of  known  sources  made 
with  which  the  standard  representation  of  the  ICRF  were  as¬ 
signed  to  the  stellar  catalog  Hipparcos  (Perryman  et  al.  1997). 
Currently  several  extensions  of  the  Hipparcos  Catalog  Reference 
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Frame  (HCRF),  densify  this  representation,  reaching  fainter 
magnitudes. 

The  quantity  of  optically  recognized  bright  quasars  has  in¬ 
creased  significantly  in  recent  years,  along  with  an  improvement 
of  the  evenness  of  their  sky  distribution  (Veron-Cetty  &  Veron 
2006;  Souchay  et  al.  2008a,b).  At  the  same  time,  the  science’s 
demands  for  accurate  astrometry  of  these  objects  have  also  in¬ 
creased.  Examples  include  micro  and  macro  lensing,  binaries, 
and  space  density  counts,  as  well  as  the  requirements  of  space 
astronomy  missions  (Andrei  et  al.  2008)  and  the  very  improve¬ 
ment  of  the  stellar  catalogs,  that  requires  for  a  dense  mesh  of 
fiducial  points  (Fienga  &  Andrei  2002,  2004). 

Therefore,  we  combine  the  large  number  of  quasars  optically 
recognized,  with  the  large,  precise  optical  catalogs  to  produce 
an  optical  materialization  of  the  ICRS  in  terms  of  its  first  princi¬ 
ples,  namely  kinematically  non-rotating  with  respect  to  the  en¬ 
semble  of  distant  extragalactic  objects,  aligned  with  the  mean 
equator  and  dynamical  equinox  of  J2000,  and  realized  by  a  list  of 
adopted  coordinates  of  extragalatic  sources.  This  article  presents 
one  such  realization,  hereafter  termed  Large  Quasar  Reference 
Frame  (LQRF). 

In  the  next  section,  the  input  data  is  presented  and  reviewed, 
whereas  in  its  final  subsection  a  schematic  of  the  data  stream 
is  presented.  The  data  stream  is  detailed  in  Sects.  3-5,  where 
the  analytical  expressions  used  to  derive  the  catalog’s  interme¬ 
diary  steps  are  shown.  In  Sect.  3,  a  local  astrometric  solution  is 
employed  to  place  all  quasar  input  catalogs  onto  the  UCAC  ref¬ 
erence  frame.  In  Sect.  4,  the  positions  derived  in  this  way  are 
aligned  with  the  ICRF  frame.  In  Sect.  5,  local  inhomogeneities 
are  tackled  by  spherical  harmonics  fitting  and  local  average  re¬ 
dressing.  The  LQRF  catalog  is  presented  in  Sect.  6,  and  a  sum¬ 
mary  and  follow-up  perspectives  are  presented  in  the  last  section. 


2.  Input  data  and  data  flow 

The  LQRF  is  built  from  gathering,  comparison,  and  adjustment 
of  positions  of  quasars  and  their  stellar  neighborhood,  extracted 
from  several  catalogs  and  lists.  A  brief  summary  of  their  content 
is  presented  in  the  following  subsections.  The  role  of  each  one 
is  also  mentioned,  for  the  sake  of  introducing  the  general  flow 
of  the  data  treatment.  The  final  subsection  outlines  the  data  flow, 
leaving  the  details  for  the  subsequent  main  sections. 

2.1.  LQAC 

The  Large  Quasar  Astrometric  catalog  (Souchay  et  al.  2008a,b), 
hereafter  LQAC,  compiles  the  detections  of  113  653  quasars 
from  the  original  entries  of  the  eight  largest  lists,  with  a  surplus 
extracted  from  the  latest  edition  of  the  Veron  &  Veron  catalog 
(Veron-Cetty  &  Veron  2006)  after  an  analysis  on  the  reliabil¬ 
ity  of  the  input  contributors.  The  LQAC  contains  extensive  in¬ 
formation  on  the  physical  properties  of  each  object.  It  contains, 
in  particular,  whenever  available,  the  redshift,  the  magnitudes 
of  the  objects  in  9  different  optical  bands  («,  b,  v,  g,  r,  i,  z,  J, 
and  K),  and  the  radio  fluxes  at  5  different  frequencies  (1.4  Ghz, 
2.3  Ghz,  5.0  Ghz,  8.4  Ghz,  and  24  Ghz).  The  position  of  each 
object  is  given  by  the  most  precise  equatorial  coordinate  directly 
available  from  the  input  lists.  Therefore,  the  precision  is  quite 
variable,  ranging  from  sub-milli-arcsec  for  VLBI  determinations 
to  a  tenth  of  arcsec  for  optical  determinations,  and  nearly  one 
arcsec  for  single  dish  determinations.  On  the  other  hand,  there 
are  no  double  entries,  either  from  misrecognition  or  from  dou¬ 
ble  radio  spots.  Although  the  position  unevenness  prevents  them 


Fig.  1.  Sky  distribution  (equatorial  coordinates,  N  up,  E  right)  of  the 
quasars  found  in  the  LQAC  catalog.  The  highest  density  regions  indi¬ 
cate  the  SDSS  DR5  contrihution.The  scale  represents  the  number  of 
sources  in  regions  of  9  sq°  shown  on  the  plots. 


Fig.  2.  Sky  distribution  (equatorial  coordinates,  N  up,  E  right)  of  the 
quasars  found  in  the  USNO  Bl.O  catalogs.  Again,  the  highest  density 
regions  indicate  the  SDSS  DR5  contrihution.The  scale  represents  the 
number  of  sources  in  regions  of  9  sq°  shown  on  the  plots.  Notice  that 
the  scale  is  slightly  different  from  that  of  Fig.  1  relative  to  the  LQAC. 


from  characterizing  an  astrometric  reference  frame,  the  LQAC 
positions  and  magnitudes  are  adopted  as  the  starting  point  for 
finding  the  objects  in  the  dense  optical  catalogs  presented  next. 
Figure  1  presents  the  LQAC  sky  distribution  and  sky  density. 

2.2.  USNO  Bl.O 

The  US  Naval  Observatory  Bl.O  catalog  (Monet  et  al.  2003), 
hereafter  Bl.O,  is  the  latest  one  of  the  series  of  all-sky  ultra 
dense  catalogs  issued  by  the  USNO.  It  brings  positions  and 
three  colors  magnitudes  (Johnston  B,  R  and  I)  at  two  epochs 
for  1  042618  261  objects,  up  to  about  V  -  20.  It  results  from 
scans  of  7435  Schmidt  plates  taken  for  six  sky  surveys  and  cali¬ 
brated  using  Tycho-2  stars  (Hog  et  al.  2000).  The  nominal  astro- 
metric  accuracy  is  200  mas,  but  the  individual  objects  position 
precision  has  a  median  value  of  1 20  mas  along  each  direction  of 
measurement.  Noticing  that  the  quasars  forming  the  LQRF  seat 
at  the  Bl.O  fainter  end,  where  astrometric  precision  and  accu¬ 
racy  is  lower,  the  LQAC  objects  were  searched  for  in  the  B 1 .0, 
and  a  neighborhood  of  stars,  i.e.,  their  position,  proper  motion, 
and  magnitude,  was  also  collected  around  the  detected  objects. 
Figure  2  presents  the  sky  distribution  and  sky  density  of  the 
83  980  quasars  recognized  in  the  USNO  Bl.O. 
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Fig.  3.  Sky  distribution  (equatorial  coordinates,  N  up,  E  right)  of  the 
quasars  found  in  the  GSC2.3  catalog.  As  before,  the  highest  density  re¬ 
gions  mark  the  SDSS  DR5  contribution.The  scale  represents  the  num¬ 
ber  of  sources  in  regions  of  9  deg^  shown  on  the  plots.  Notice  that  the 
scale  is  slightly  different  from  that  of  Fig.  2  relative  to  the  B 1 .0  catalog. 


2.3.  GSC2.3 

The  Second  Generation  Guide  Star  catalog  latest  version  (Lasker 
et  al.  2008),  hereafter  GSC2.3,  is  an  all-sky  catalog  derived  from 
the  Digitized  Sky  Survey  (Taff  et  al.  1990a,b)  that  the  Space 
Telescope  Science  Institute  and  the  Osservatorio  di  Torino  cre¬ 
ated  from  the  Palomar  and  UK  Schimdt  survey  plates.  It  con¬ 
tains  position,  proper  motion,  and  magnitude  (Johnston  B,  R, 
and  I)  for  945  592  853  objects,  and  is  expected  to  be  complete 
to  R  =  20.  The  astrometric  total  error  is  quoted  as  smaller 
than  300  mas,  while  the  relative  astrometric  error  is  lower  than 
200  mas.  In  fact,  the  astrometric  residuals  compared  to  exter¬ 
nal  reference  catalogs  only  become  higher  than  150  mas  above 
R  -  1 8  at  the  fainter  end.  The  astrometric  calibration  was  done 
using  Tycho-2  stars.  As  can  be  seen,  there  are  overwhelming 
similarities  between  the  GSC2.3  and  the  Bl.O  raw  data,  which 
were  measured  and  treated  however  by  distinct  methods,  en¬ 
abling  a  most  natural  combination  of  the  independent  results 
from  each  of  them.  As  before,  we  searched  for  the  LQAC  objects 
in  the  GSC2.3,  and  a  neighborhood  of  stars  (position,  proper  mo¬ 
tion,  and  magnitude)  was  collected  around  the  found  objects. 
Figure  3  presents  the  sky  distribution  and  sky  density  of  the 
93  943  quasars  recognized  in  the  GSC2.3. 

2.4.  SDSS  DR5 

The  Sloan  Digitized  Sky  Survey  (SDSS)  data  release  5 
(Adelman-McCarthy  et  al.  2007),  hereafter  DR5,  represents  the 
completion  of  the  SDSS-I  project.  It  covers  8000  sq°,  around 
the  north  Galactic  cap  (b  <  30°)  and  three  stripes  towards  the 
south  Galactic  cap,  for  magnitudes  as  faint  as  i  =  20.2,  and  a 
total  of  215  million  unique  objects.  The  astrometry  is  prefer¬ 
ably  referred  to  the  UCAC  catalog,  or  else  to  the  Tycho-2  cat¬ 
alog.  The  rms  residuals  per  coordinate  are  45  mas  and  75  mas 
for  reductions  against  the  UCAC  (Zacharias  et  al.  2000)  and  the 
Tycho-2  (Hoeg  et  al.  1997),  respectively,  with  additional  30  mas 
systematic  errors  in  both  cases  (Pier  et  al.  2003).  The  internal 
astrometric  precision  is  mainly  limited  by  seeing  effects,  and  it 
is  quoted  as  around  35  mas  in  each  coordinate.  We  note  that 
quasars  are  signaled  in  the  DR5,  for  a  total  of  74  869  entries 
with  quasi-stellar  object  spectral  class.  As  seen  from  Fig.  1,  the 
large  number  of  DR5  quasars  represents  an  important  feature  of 
the  LQAC  input. 


2.5.  UCAC2 

The  second  USNO  CCD  Astrograph  catalog  (Zacharias  et  al. 
2004),  hereafter  UCAC2,  is  the  latest  release  of  the  ongoing 
UCAC  project,  designed  to  observe  the  entire  sky  for  R  mag¬ 
nitudes  between  about  7.5  and  16.  The  observed  positional  er¬ 
rors  are  about  20  mas  for  the  stars  in  the  10  to  14  mag  range, 
and  about  70  mas  at  the  limiting  magnitude  of  R  -  16.  The 
UCAC2  is  a  high  density,  highly  accurate,  astrometric  catalog 
(positions  and  proper  motions)  of  48  330571  stars  covering  the 
sky  from  -90  to  -1-40°  in  declination  and  as  high  as  h-52°  in  some 
areas.  The  northern  limit  is  a  function  of  right  ascension.  The 
astrometry  provided  in  the  UCAC2  is  on  the  Hipparcos  system. 
Positions  are  given  at  the  standard  epoch  of  Julian  date  2000.0. 
The  UCAC2  supersedes  the  Tycho-2  catalog,  having  a  stellar 
density  that  is  over  20  times  higher.  It  is  used  here  as  the  pre¬ 
ferred  catalog  to  place  the  quasar  positions  collected  from  the 
catalogs  described  in  the  above  subsections  onto  the  Hipparcos 
system.  Nevertheless,  since  the  UCAC2  coverage  does  not  in¬ 
clude  the  entire  sky,  two  additional  stellar  frames  are  comple¬ 
mentary  used  and  are  presented  below. 

2.6.  UCACN 

We  hereafter  term  UCACN  a  cutout  of  the  preliminary  UCAC 
positions  around  the  LQAC  quasars  in  the  northern  part  of  the 
sky  that  is  not  reached  by  the  UCAC2  coverage  (Zacharias 
2007).  Although,  as  discussed  further  on,  the  astrometric  accu¬ 
racy  of  the  UCACN  positions  and  its  stellar  density  granted  no 
degradation  on  the  results  obtained  using  it,  it  will  be  dropped 
in  future  revisions  of  this  work,  after  the  forthcoming  release  of 
the  all-sky  catalog  UCAC3.  The  nominal  astrometric  accuracy 
of  the  UCACN  stars  is  similar  to  that  of  the  UCAC2  ones,  but 
no  proper  motion  data  were  released  for  these  stars.  Throughout 
this  paper  when  referring  to  the  combination  of  the  UCAC2  and 
UCACN,  the  designation  UCAC  is  used. 

2.7.  2MASS 

The  Two  Micron  All-Sky  Survey  point  source  catalog  (Cutri 
et  al.  2003),  hereafter  2MASS,  derives  from  an  uniform  scan 
of  the  entire  sky  in  three  near-infrared  bands  to  detect  and  char¬ 
acterize  point  sources  brighter  than  about  1  mJy  in  each  band, 
with  signal-to-noise  ratio  greater  than  10,  using  a  pixel  size 
of  2.0",  each  point  in  the  survey  having  been  imaged  six  times. 
The  detectors  worked  to  a  3cr  limiting  sensitivity  of  J  -  17.1, 
H  -  16.4,  and  K  -  15.3.  The  2MASS  contains  the  position  of 
470992970  sources,  but  no  proper  motions.  The  astrometry  is 
referred  to  the  Tycho-2  catalog  and  it  is  accurate  to  70-80  mas 
over  the  magnitude  range  of  9  <  K  <  14  mag.  Comparative 
studies  have  shown  that  the  2MASS  positions  agree  with  those 
from  common  UCAC2  stars  to  within  10  mas  (Zacharias  et  al. 
2003),  and  that  the  2MASS  is  compliant  to  the  ICRF  in  the  range 
100-120  mas.  Taking  advantage  of  its  all-sky  coverage,  very 
high  stellar  density,  and  consistent  astrometry,  the  2MASS  cat¬ 
alog  is  used  here  as  an  independent  means  to  place  the  LQAC 
quasars  recognized  in  the  large,  deep  catalogs  on  the  Hipparcos 
system. 

2.8.  ICRF-Ext2 

The  International  Celestial  Reference  Frame  2nd  extension  (Ma 
et  al.  1998;  Fey  et  al.  2004),  hereafter  ICRF-Ext2,  is  the  present 
materialization  of  the  International  Celestial  Reference  System 


388 


A.  H.  Andrei  et  al.:  The  LQRF  -  an  optical  representation  of  the  ICRS 


(ICRS)  at  radio  frequencies.  It  represents  the  basic  frame  with 
respect  to  which  the  position  of  any  object  in  the  celestial  sphere 
should  be  measured.  In  its  primary  contents,  the  ICRF  con¬ 
sists  of  212  sources,  called  “defining”,  whose  positions  are  inde¬ 
pendent  of  the  classical  planes  (equator,  ecliptic)  and  reference 
points  (equinox),  but  consistent  with  the  previous  realizations 
of  the  Celestial  System  as  the  FK5.  With  the  help  of  a  con¬ 
siderable  amount  of  VLBI  observations,  the  individual  posi¬ 
tions  of  the  sources  were  found  to  be  accurate  to  within  roughly 
0.25  mas,  while  the  stability  of  the  reference  axes  attains  a  re¬ 
markable  20  fias  (micro-arcsec)  accuracy.  In  the  ICRF  initial 
version  608  sources  complemented  the  defining  212  sources. 
The  entrance  of  the  ICRF-Extl  and  of  the  ICRF-Ext2  provided 
109  new  sources  to  the  catalog,  thus  leading  to  a  total  number 
of  717  radio-sources.  Notice  that  in  the  ICRE-Ext.2  catalog,  a 
very  small  sample  of  objects  are  not  quasars:  10  of  them  are 
cathegorized  as  AGN  (Active  Galactic  Nuclei)  and  10  are  cathe- 
gorized  as  BE  LAC  (BE  Lacertae).  Although  the  LQAC  a  priori 
exclusively  considers  quasars  in  its  compilation,  these  particular 
objects  were  retained  because  of  their  astrometric  accuracy.  The 
ICRE-Ext2  is  used  here  to  match  the  VLBI  positions  directly  to 
the  derived  optical  positions.  In  the  derivation  of  the  LQRE,  such 
match  serves  to  two  distinct  steps:  (a)  to  redress  the  global  frame 
orientation  towards  the  ICRF;  and  (b)  to  redress  local  departures 
of  the  derived  quasar  catalog  from  the  ICRF. 


2.9.  VCS6 

The  Very  Long  Baseline  Array  (VLBA)  Calibrator  Survey 
(Petrov  et  al.  2006),  hereafter  VCS6,  consists  of  a  catalog  con¬ 
taining  milli-arcsec  accurate  positions  of  3910  extragalacic  radio 
sources,  mainly  quasars.  These  positions  were  derived  from  the 
astrometric  analysis  of  dual-frequency  2.3  and  8.4  GHz  VLBA 
observations,  in  the  framework  of  the  International  Celestial 
Reference  System  (ICRS)  as  realized  by  the  ICRF.  Thus,  to 
the  level  of  precision  and  accuracy  of  the  optical  positions  de¬ 
termined  here,  the  VCS6  sources  represent  the  ICRF.  In  order 
to  strictly  validate  such  representation,  524  VCS6  sources  with 
formal  inflated  errors  larger  than  10  mas  were  excluded.  The 
VCS6  sources  that  did  not  belong  to  the  ICRF-Ext2  were  added 
to  it,  to  obtain  an  enlarged  radio  frame,  which  was  used  for  the 
two  steps  mentioned  in  the  ICRF-Ext2  subsection,  that  is,  the 
global  and  local  redressing  of  the  optical  quasar  positions  to¬ 
wards  the  ICRF. 


2.10.  VLAC 

The  Very  Large  Array  (VLA)  Calibrator  list  (Claussen  2006), 
hereafter  VLAC,  consists  of  5523  radio  flux,  structure  code, 
position,  and  associated  error  entries  for  1860  extragalactic 
sources.  It  is  foremostly  indicated  as  a  reference  tool  for  ob¬ 
servations  at  the  VLA,  but  the  radio  positions  are  accurate  and 
compliant  with  the  ICRF.  Then,  as  in  the  previous  subsection, 
the  VLAC  positions  were  added  to  those  of  the  ICRF-Ext2  core 
data  set,  after  removing  297  with  errors  larger  than  10  mas.  The 
combination  of  the  ICRF-Ext2,  VCS6,  and  VLAC  will  be  here 
termed  as  the  enlarged  radio  frame  (ERF)  and  represents  the 
ICRF  to  within  an  uncertainty  of  10  mas.  The  order  of  prior¬ 
ity  for  the  radio  lists  contributing  to  the  enlarged  radio  frame 
was,  ICRF-Ext2  (718  sources),  VCS6  (2684  sources),  and  Anally 
VLAC  (59  sources).  The  3461  sources  now  represent  the  ICRF 
on  a  better  distributed  all-sky  basis  (Fig.  4),  and  the  sources  for 


Fig.  4.  Sky  distribution  (equatorial  coordinates,  N  up,  E  right)  of  the 
enlarged  radio  frame  quasars.  The  filled  symbols  mark  the  sources  for 
which  an  optical  counterpart  is  found  in  at  least  one  of  the  deep,  dense 
source  catalogs  (Bl.O,  GSC2.3,  DR5).  The  enlarged  radio  frame  is 
formed,  in  order  of  priority  of  choice,  by  the  ICRF-Ext2  (718  sources), 
the  VCS6  (2684  sources),  and  the  VLAC  (59  sources). 


which  the  optical  counterpart  position  is  known  will  provide  the 
most  reliable  means  of  relating  the  LQRF  to  the  ICRF. 

2.11.  Data  flow 

The  various  origins  of  the  sources  compiled  in  the  LQAC  im¬ 
poses  first  of  all  to  refer  their  positions  to  a  single  frame  on  the 
ICRS  system.  After  that,  the  quasar  frame  formed  in  this  way 
is  aligned  with  the  ICRF.  Finally,  the  position  of  each  source  is 
calculated  by  weighted  average.  This  data  flow  can  be  sketched 
as  follows: 

1 .  The  LQAC  entries  are  all  admitted. 

2.  The  LQAC  quasars  are  searched  for,  by  a  main  criterion  of 
position  coincidence,  in  the  Bl.O,  GSC2.3,  and  SDSS  DR5 
high  density  catalogs.  Each  individual  match  is  kept  and 
treated  separately. 

3.  Around  each  of  the  individual  quasars  thus  matched  a  stellar 
neighborhood  is  extracted  from  the  appropriate  high  density 
catalogs. 

4.  The  stellar  neighborhood  contents  from  the  high  density 
catalogs  are  searched  for  in  the  stellar  catalogs  UCAC2, 
UCACN,  and  2MASS,  as  before  by  a  main  criterion  of  posi¬ 
tion  coincidence. 

5.  Within  each  neighborhood,  the  stars  common  to  the  high 
density  catalog  and  to  the  stellar  catalog  have  their  positions 
projected  onto  the  tangential  plane  centered  on  the  quasar.  A 
local  transformation  function  from  the  high  density  catalog 
to  the  stellar  catalog  is  calculated,  mimicking  the  complete 
plate  solution  polynomials. 

6.  At  this  point,  for  each  quasar,  typically,  a  number  of  posi¬ 
tions  are  determined  (e.g.,  from  more  than  one  high  den¬ 
sity  catalog,  and  local  solutions  from  the  two  stellar  cata¬ 
logs).  Each  of  these  families  of  positions  will  still  be  handled 
separately. 

7.  Within  each  family,  the  quasars  for  which  there  is  a  position 
in  the  enlarged  radio  frame  (that  is,  the  combination  of  the 
ICRE,  VCS6,  and  VLAC  catalogs)  are  found,  since  they  are 
also  flagged  in  the  LQAC,  what  as  before  follows  mainly  a 
criterion  of  positional  coincidence. 
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8.  The  subsets  so  defined  are  used  to  calculate  the  global  ro¬ 
tation  and  zero  point  corrections  to  obtain  coincidence  with 
the  ICRF  origin  and  axes.  The  global  rotation  and  zero  point 
corrections  are  next  applied  to  the  positions  of  all  quasars  of 
each  given  family. 

9.  As  in  the  previous  step,  the  optical  and  radio  positions  in  sub¬ 
sets  are  combined  using  orthogonal  functions,  in  right  ascen¬ 
sion,  declination,  and  magnitude,  to  correct  systematic  local 
departures  from  the  ICRF  axes  directions.  And  likewise  the 
significant  orthogonal  functions  are  applied  to  the  positions 
of  all  quasars  in  each  given  family. 

1 0.  Again  the  optical  and  radio  positions  in  subsets  are  combined 
in  small  regions  around  each  quasar  to  derive  its  correction 
for  localized  inhomogeneities. 

1 1 .  For  each  position  of  each  source  in  a  given  family,  a  total 
error  is  assigned  that  combines  the  formal  error  from  the  high 
density  catalog  where  it  appears,  the  error  from  the  particular 
local  correction,  the  error  from  the  particular  global  rotation 
and  bias,  and  the  error  from  both  the  significant  orthogonal 
function  adjustments  and  the  inhomogeneity  corrections. 

12.  From  the  squared  inverse  of  the  total  error,  a  weight  is  as¬ 
signed  to  each  determined  position  of  a  given  quasar.  The 
final  position  of  each  quasar  is  given  as  the  weighted  aver¬ 
age  of  the  determined  positions.  The  error  in  the  position  is 
assigned  according  to  an  individualized  function  of  the  opti¬ 
cal  to  radio  disagreement. 

The  comparisons  of  the  final  positions  against  the  existing  cor¬ 
responding  positions  in  the  enlarged  frame  is  used  to  obtain  the 
global  properties  of  the  LQRF. 


3.  The  local  astrometric  solution 

We  now  derive  the  points  1  to  6  of  Sect.  2.1 1  (data  flow). 

Starting  from  the  LQAC  input  list  of  quasars,  the  first  step 
in  obtaining  their  astrometric  positions,  to  form  a  whole  consis¬ 
tent  reference  frame,  is  to  match  the  positions  of  the  objects  with 
those  in  the  two  large,  deep  all  sky  catalogs,  the  USNO  B 1 .0  and 
the  GSC2.3,  as  well  as  in  the  SDSS  DR5.  The  basic  matching 
criteria  is  position  agreement  to  within  1  arcsec.  This  is  equal 
to  the  poorest  astrometric  accuracy  in  the  surveys  contribut¬ 
ing  to  the  LQAC.  It  is  also  a  safe  threshold  given  the  typical 
seeing  and  astrographic  plate  resolution  limits.  Additionally,  a 
magnitude  limit  was  applied  preventing  the  selection  of  objects 
brighter  than  R  -  1 ,  nevertheless  this  threshold  was  not  reached 
by  any  source.  No  proper  motion  limit  was  enforced  because 
both  the  USNO  Bl.O  and  the  GSC2.3  warn  against  possible  off¬ 
sets  of  the  proper  motion  zeropoints.  A  check  was  made  after¬ 
wards  by  determining  local  zeropoints  of  proper  motion,  given 
by  the  average  motion.  In  comparison  to  the  local  zeropoints  of 
proper  motion,  just  0.7%  of  the  GSC2.3  selected  objects  show 
significant  motion,  while  for  the  Bl.O  selected  objects  just  0.1% 
of  them  show  significant  motion.  To  verify  what  would  be  the 
odds  of  a  chance  hit,  the  input  coordinates  were  varied  by  ran¬ 
dom  values  between  1  and  5  arcmin.  In  this  way  the  same  re¬ 
gion  of  the  true  objects  was  swept.  Keeping  the  search  radius 
to  1  arcsec,  in  only  0.2%  of  the  cases  was  verified  a  hit  on 
the  false  coordinates.  Moreover,  in  only  0.1%  of  cases  was  a 
second  object  found  significantly  close  to  the  adopted  match. 
These  satisfying  rates  change  very  little  when  the  search  radius 
is  increased  to  2  arcsec.  While  this  provides  additional  proof 
of  the  overwhelming  correctness  of  the  matches,  it  at  the  same 
time  indicates  that  there  is  little  advantage  in  adopting  a  larger 


Table  1.  Initial  optical  to  radio  relationships" . 


Tracer 

Bl.O  to  RadPos 

GSC2.3  to  RadPos 

DR5  to  RadPos 

N 

1962 

2028 

300 

Aa  cos  S 

-35.2  ±  4.7 

+32.9  ±5.5 

^.5  ±  3.2 

A6 

+  123.5  ±4.9 

+39.5  ±  4.7 

+  13.9  ±3.2 

Al 

+4.7  ±  6.2 

-6.6  ±  6.7 

-3.5  ±  4.7 

A2 

+5.9  ±5.8 

+  13.1  ±6.3 

-19.4  ±5.1 

A3 

+42.1  ±  5.7 

-29.0  ±6.1 

+6.6  ±  4.1 

A4N 

+  123.0  ±5.3 

20.9  ±  5.7 

1.6  ±4.5 

A4S 

+  128.0  ±  11.9 

134.1  ±  12.5 

“  For  the  catalogs  USNO  Bl.O,  GSC2.3,  and  SDSS  DR5  relatively 
to  the  collected  sample  of  sources  with  precise  radio  interferometric 
positions  (Extended  Radio  Frame  -  RadPos,  either  from  ICRF-Ext2, 
VCS6  or  VLAC),  are  presented:  the  average  offsets  in  right  ascension 
and  declination,  the  orientation  angles  (A1-A3)  relative  to  the  standard 
thriedron  of  equatorial  coordinates,  and  the  equatorial  biases  (A4N  and 
A4S,  north  and  south  of  -20°).  All  values  are  in  milli-arcsec  (mas)  and 
the  number  of  sources  is  given  in  the  first  line. 


radius,  and  certainly  no  gain  to  the  confidence  level  attached  to 
the  cross-identifications. 

In  all,  83  980  quasars  are  recognized  in  the  Bl.O,  and 
93  943  quasars  are  recognized  in  the  GSC2.3.  The  sky  distri¬ 
bution  follows  closely  that  shown  in  Fig.  1 .  From  the  figure  the 
significant  contribution  of  the  SDSS  DR5  is  clear.  For  the  DR5, 
74  825  quasars  are  retrieved  using  the  SDSS  query  server.  The 
positions  from  the  three  catalogs  were  added  independently  to 
the  input  list.  The  final  input  data  contains  100  165  quasars.  Just 
small  fractions  of  the  total  appear  in  only  one  of  the  catalogs 
(1405  only  in  USNO  Bl.O,  4425  only  in  GSC2.3,  and  3276  only 
in  SDSS  DR5),  while  61  532  quasars  appear  in  all  three  catalogs. 

These  three  catalogs  that  bring  the  input  positions  for  the 
recognized  quasars  are  originally  placed  on  the  ICRF  J2000 
reference  frame  by  different  astrometric  pipelines  based  on  the 
Tycho-2  catalog.  Their  global  orientation  and  equatorial  bias  rel¬ 
ative  to  the  ICRF  J2000  frame  are  shown  in  Table  1.  Local  devi¬ 
ations  are  also  apparent  in  Figs.  5-7. 

To  homogenize  all  the  original  different  positions,  a  local  as¬ 
trometric  reduction  was  performed  based  on  the  catalogs  UCAC 
and  2MASS.  The  local  correction  follows  the  standard  tangent 
plane  astrometric  reduction,  taking  the  quasar  original  position 
as  the  central  point  of  the  rigorous  gnomic  projection  (Assafin 
et  al.  1997).  Through  it  the  catalog  coordinates  are  adjusted  to 
the  reference  frames  represented  by  the  stellar  catalogs.  Three 
assumptions  are  implicit  to  the  method:  that  the  stellar  catalogs 
deliver  an  improved,  and  similar,  representation  of  the  ICRF,  ei¬ 
ther  relatively  to  the  two  Schmidt  plate  surveys,  as  well  as  rel¬ 
atively  to  the  photometric  one.  That  in  a  very  restricted  neigh¬ 
borhood  the  deviation  between  the  deep  surveys  and  the  stellar 
catalogs,  whichever  one  it  originated  from,  can  be  modeled  by 
simple  relationships.  And  that  the  removal  of  such  deviations 
can  significantly  improve  the  accuracy  of  the  quasar  positions 
found  in  the  large  surveys.  The  results  of  previous  investigations 
support  these  assumptions  (da  Silva  Neto  et  al.  2005). 

3.1.  The  polynomial  degree  of  the  solution 

On  the  functional  side,  the  initial  question  to  be  tackled  is  the 
polynomial  degree  of  the  local  astrometric  correction  represen¬ 
tation.  The  UCAC  is  the  less  dense  of  the  catalogs  involved,  and 
already  presents  a  typical  density  of  10  stars  per  5x5  arcmin. 
Translated  to  the  Schmidt  plates  scale  this  means  that  a  first 
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Fig.  5.  USNO  Bl.O  right  ascension  and  declination  offset  distributions  relative  to  the  equatorial  coordinates  axes  and  the  R  magnitude.  The  offsets 
are  given  in  the  sense  of  catalog  minus  the  radio  interferometric  positions  from  the  ICRF-Ext2,  or  the  VCS6,  or  the  VLAC. 
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Fig.  6.  GSC2.3  right  ascension  and  declination  offsets  distributions  relative  to  the  equatorial  coordinates  axes  and  the  R  magnitude.  The  offsets  are 
given  in  the  sense  of  catalog  minus  the  radio  interferometric  positions  from  the  ICRF-Ext2,  or  the  VCS6,  or  the  VLAC. 
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Fig.  7.  SDSS  DR5  right  ascension  and  declination  offsets  distributions  relative  to  the  equatorial  coordinates  axes  and  the  R  magnitude.  The  offsets 
are  given  in  the  sense  of  catalog  minus  the  radio  interferometric  positions  from  the  ICRF-Ext2,  or  the  VCS6,  or  the  VLAC. 


degree  solution  with  10  reference  stars  can  be  applied  over  a 
physical  area  as  small  as  25  mm^.  According  to  the  adopted  a 
priori  assumptions,  and  given  the  sky  density  of  the  stellar  cata¬ 
logs,  the  simplest  models  seem  to  be  well  suited.  To  verify  this 
procedure,  a  trial  local  reduction  was  performed  for  29  quasars 
randomly  distributed  across  the  celestial  sphere,  at  every  one 
hour  of  right  ascension,  and  away  from  the  galaxy  plane.  Around 
each  of  them,  a  field  of  UCAC2  stars  was  selected  in  boxes 
of  side  10  arcmin,  which  contained  at  least  30  stars.  For  these 
29  fields  of  quasars,  complete  polynomial  solutions  were  tested, 
with  orders  varying  from  0  to  3rd,  and  a  straight  averages  solu¬ 
tion.  The  quasar  positions  were  taken  from  the  USNO  Bl.O  and 
all  positions  were  transported  to  the  USNO  Bl.O  plate  epoch, 
using  the  proper  motions  of  either  catalog.  Table  2  summarizes 
the  mean  results. 

In  the  hrst  half  of  Table  2,  the  average  local  corrections  and 
the  corresponding  standard  deviation  are  shown.  Although  the 
values  only  very  loosely  represent  the  systematic  offsets  between 
the  USNO  Bl.O  and  the  UCAC2,  it  is  important  to  note  is  the 
similarity  between  all  the  rows.  In  the  second  half  of  the  ta¬ 
ble,  the  average  internal  errors  of  each  solution  and  their  dis¬ 
persion  are  shown.  Again  all  rows  behave  quite  alike.  Therefore, 
to  minimize  the  region  in  which  the  local  solution  is  derived, 
and  maximize  the  ratio  of  the  number  of  equations  to  the  num¬ 
ber  of  parameters  to  find,  the  complete  hrst  degree  polynomial 
was  adopted.  This  solution  also  produces  the  smallest  set  of  stan¬ 
dard  deviations.  The  zero  degree  solutions  fared  just  as  well  as 
the  others,  so  there  would  be  no  particular  advantage  in  adopting 
an  even  smaller  region,  while  the  hrst  degree  polynomial  allows 
for  some  modeling,  which  is  useful  face  to  the  large  numbers  of 


Table  2.  Comparison  between  the  models  for  the  local  astrometric 
solution". 


Solution 

Aa  cos  S 

A6 

O'  AS 

3rd  degree 

+A2 

126 

-25 

140 

2nd  degree 

+43 

127 

-23 

136 

1st  degree 

+61 

110 

-34 

120 

0  degree 

+42 

126 

-34 

119 

offset 

+56 

118 

-45 

130 

Solution 

O'l.x 

Cey 

3rd  degree 

137 

63 

127 

64 

2nd  degree 

134 

66 

125 

64 

1st  degree 

135 

68 

120 

65 

0  degree 

133 

68 

117 

64 

“  In  the  top  half  of  the  table,  the  average  corrections  and  their  stan¬ 
dard  deviations  are  presented.  In  the  bottom  half,  the  average  internal 
error  in  the  solutions  and  their  dispersions  are  presented.  All  values  are 
in  milli-arcsecs.  The  difference  between  the  “0  degree”  and  the  “off¬ 
set”  solutions  is  just  that  in  the  former  residuals  larger  than  2.5cr  are 
removed  and  the  solution  is  recalculated. 


sources  that  the  local  solution  is  going  to  be  applied  to.  Two  for¬ 
mulations  of  the  hrst  degree  polynomial  are  concurrently  used, 
the  complete  one  (six  independent  parameters)  -  hereafter  the 
six-parameter  model,  and  the  so-called  four  constants  one  (in 
which  the  X  and  Y  axes  are  constrained  to  have  the  same  scale), 
hereafter  the  four-parameter  model.  Furthermore,  we  note  that 
with  the  six-parameter  model  the  error  in  the  central,  pseudo- 
tangential  point,  which  is  always  the  position  of  the  quasar,  was 
also  found  to  be  the  smallest. 
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3.2.  The  stellar  catalogs  and  the  size  of  the  local  segion 

As  discussed  above,  the  UCAC2  provides  the  desiderata  for  the 
local  astrometric  solution:  high  stellar  density,  accurate  global 
representation  of  the  HCRF  and  hence  of  the  ICRF,  range  of 
magnitudes  up  to  16,  and  small  zonal  and  local  frame  biases. 
The  stellar  density  supports  using  a  region  as  small  as  10  arcmin 
across  to  be  used,  which  is  always  within  the  overlap  zone  of  the 
Schmidt  plate  surveys  contributing  to  the  Bl.O  and  the  GSC2.3. 

However,  the  UCAC2  has  incomplete  coverage  in  the  north¬ 
ern  celestial  hemisphere.  Two  remedies  are  used  to  account 
for  this  hindrance,  the  use  of  preliminary  positions  of  northern 
UCAC  stars  (here  called  UCACN),  and  the  use  of  the  2MASS 
catalog.  In  the  former  case,  the  UCACN  positions  are  by  con¬ 
struction  on  the  same  system  of  the  UCAC2,  and  of  the  same 
precision. 

The  2MASS  Point  Source  catalog  positions  are  compliant 
with  the  ICRS  by  means  of  the  Tycho-2  catalog,  and  are  expected 
to  be  accurate  to  70-80  mas  over  the  magnitude  range  of  9  < 
Ks  <  14  mag.  At  the  limit  Ks  =  16,  the  accuracy  is  at  the  level  of 
200  mas.  As  the  2MASS  does  not  provide  proper  motions,  these 
values  are  good  for  the  mean  epoch  of  observation,  1999. 

The  pitfalls  created  by  the  magnitude  gap  between  the 
2MASS  objects  and  the  brighter  Tycho-2  reference  stars  re¬ 
duce  the  final  astrometric  accuracy.  The  internal  precision  of  the 
2MASS  individual  positions  is  better,  at  the  level  of  40-50  mas, 
in  the  9  <  ^fs  <  14  mag  range.  A  comparison  with  5.2  million 
stars  also  present  in  the  UCAC2  shows  standard  deviation  of 
70-80  mas  in  both  coordinates,  for  the  9  <  TTs  <  14  mag  range. 
For  the  stars  in  common  around  the  quasars  extracted  to  de¬ 
fine  the  local  solution,  the  difference  average  and  standard  de¬ 
viation  are  Aacosd  =  -i-3  mas,  cr^a  cos  s  =  99  mas,  and  Ad  = 
-Hll  mas,  (Tas  -  95  mas,  all  in  the  sense  2MASS. minus  UCAC2. 
Furthermore,  from  the  limited  sample  of  700  quasars  from  the 
enlarged  radio  frame  found  in  the  2MASS  catalog,  the  analo¬ 
gous  statistics  are  Aacosd  =  H-0  mas,  crAacosS  =  144  mas,  and 
Ad  =  -2  mas,  ctas  =  138  mas. 

On  the  positive  side,  for  the  same  region,  the  2MASS  cat¬ 
alog  contains  about  10  times  the  number  of  reference  stars  as 
the  UCAC  catalog.  In  the  small  regions  used  for  the  local  solu¬ 
tion,  it  represents  a  coherent  local  frame,  from  which  the  zonal 
bias  can  be  removed  after  later  comparisons  with  the  enlarged 
radio  frame  sources.  It  is  well  aligned  with  both  the  UCAC  and 
the  ICRF,  and  represents  an  all-sky  frame.  The  local  solution 
was  computed  independently  using  UCAC  and  2MASS  refer¬ 
ence  stars.  Each  of  the  results  contribute,  again  independently, 
to  the  weighted  average  used  to  compute  the  final  LQRF  quasar 
coordinates. 

3.3.  The  analytical  expression  for  the  local  correction 
astrometric  solution 

On  establishing  the  analytical  form  of  the  local  astrometric  so¬ 
lution,  a  final  point  to  be  considered  concerns  the  usefulness  of 
a  magnitude  dependent  term,  a  feature  often  present  in  Schmidt 
plate  based  surveys  (da  Silva  Neto  et  al.  2000).  In  the  present 
case,  the  R  magnitude  represents  the  obvious  choice,  mainly 
because  it  is  common  to  all  the  involved  catalogs,  and  closer 
to  the  central  wavelengths  of  the  2MASS  bands.  At  the  same 
time,  because  of  the  quasars  redshift  and  to  the  reddening  for  the 
stars,  coupled  with  CCD  top  sensitivity  and  wider  plate  effec¬ 
tive  passband,  the  signal-to-noise  ratio  is  usually  higher  towards 
the  red.  An  inspection  on  the  magnitude  dependency  of  the 


catalog  minus  enlarged  radio  frame  residuals  indicates  a  small 
linear  trend,  which  has  a  higher  value  for  the  A6  offsets  than 
for  the  Aa cos  6  offsets.  The  same  trend  is  verified  for  both  Bl.O 
and  CSC2,3  increasing  with  magnitude  for  declination,  and  de¬ 
creasing  with  magnitude  for  right  ascension.  The  amplitudes  are 
not  higher  than  10  mas  in  the  significant  cases.  For  the  DR5, 
the  magnitude  dependence  is  of  much  smaller  amplitude  and 
opposite  trend.  The  possibility  of  a  magnitude  equation  is  not 
excluded  in  the  presentation  of  either  of  the  two  deep  catalogs. 
However,  the  causes  and  ensuing  description  would  be  out  of 
the  scope  of  this  work.  Here,  an  additional  obstacle  would  be  to 
bridge  the  magnitude  gap  between  the  reference  stars  and  the  far 
dimmer  quasars.  All  accounted,  a  magnitude  dependency  is  left 
to  be  investigated  by  the  orthogonal  functions  in  direct  compar¬ 
ison  with  the  enlarged  radio  frame  positions. 

Therefore,  the  enforced  local  solution  is 

^SCat  -  TfQCat  =  ri  TCgCat  +  b  TgCat  +  C,  (1) 

?7SCat  -  Tgcat  =  d  iQCat  +  £  TCgcat  +  /•  (2) 

In  Eqs.  (1)  and  (2),  “SCat”  is  either  the  UCAC2  (or  even  its 
preliminary  northernmost  part)  or  the  2MASS,  and  “QCat”  des¬ 
ignates  the  Bl.O,  the  CSC2.3  or  the  DR5.  The  relations  above 
were  determined  for  the  smallest  set  of  “SCat”  and  “QCat”  stel¬ 
lar  objects  in  commonin  a  neighborhood  larger  than  10  arcmin 
but  smaller  than  a  box  of  size  of  30  arcmin,  to  a  minimum  of 
6  stars.  The  solutions  were  calculated  independently  for  right 
ascension  and  declination.  Whenever  that  minimum  number  of 
stars  is  touched,  either  at  the  start  or  after  elimination  of  stars, 
the  four-parameter  model  is  used,  combining  the  two  equato¬ 
rial  coordinate  solutions.  In  fact,  as  shown  in  the  next  section, 
the  four-parameter  model  was  applied  for  all  quasars  and  it  de¬ 
livers  solutions  in  excellent  agreement  with  the  six-parameter 
solutions.  The  analytical  form  of  the  four-parameter  model  is 
presented  below  in  Eqs.  (3)  and  (4),  where  “SCat”  and  “QCAT” 
keep  their  afore-defined  meaning. 

^SCat  -  TfQCat  =  g  A'gcat  +  h  TgCat  +  U  (3) 

'iSCat  -  Tgcat  =  h  TgCat  “  9  A'QCat  +  j-  (4) 


3.4.  Applying  the  local  astrometric  solution 

As  discussed  in  the  previous  sections,  the  initial  task  of  produc¬ 
ing  the  LQRF  consists  of  carrying  out  a  local  astrometric  re¬ 
duction  using  the  most  precise  and  accurate  stellar  catalogs  rep¬ 
resenting  the  HCRF,  through  the  Tycho-2  reference  frame,  to 
locally  redress  the  quasar  position  input  from  the  high  density, 
deep  catalogs. 

There,  however,  the  question  of  epochs  comes  up.  The  stel¬ 
lar  catalogs  have  a  mean  epoch  very  close  to  J2000,  being 
J1998.79  (+2.78  y)  for  the  UCAC2,  and  12003.34  (+0.24  y) 
for  the  UCACN,  the  UCAC2  preliminary  northernmost  part,  and 
J1998.87  (+0.08  y)  for  the  2MASS.  Also,  the  DR5  mean  epoch 
is  J2002.96  (+1.45  y),  thus  raising  no  questions.  The  dilemma  on 
whether  the  local  reduction  ought  to  be  performed  in  the  epoch 
of  the  quasar  input  position  or  directly  in  J2000  is  only  mean¬ 
ingful  when  the  Bl.O  and  the  CSC2.3  are  locally  reduced  by  the 
UCAC2.  Then,  it  would  not  be  unreasonable  to  speculate  that  a 
reduction  made  at  the  quasar  epoch  would  diminish  the  scatter 
in  the  stellar  positions,  on  the  supposition  that  the  stellar  im¬ 
ages  themselves  were  at  the  quasar  epoch  (i.e.,  the  epoch  of  the 
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Table  3.  Mean  results  from  the  local  astrometric  corrections". 


Solution 

Nq 

Ns 

Ca  cos^ 

CS 

Nrad 

Att  cos  6 

AS 

BU6 

47  822 

13.0 

-23.5  (106.2) 

-145.3  (104.2) 

1320 

^3.5  (177.0) 

+  12.2(171.6) 

BN6 

3023 

14.1 

+4.2  (45.6) 

-230.0  (81.7) 

123 

-70.5  (136.6) 

-9.9(121.6) 

BT6 

83  951 

52.8 

+3.3  (107.9) 

-178.8  (111.2) 

1897 

-18.2  (172.9) 

-12.0  (167.5) 

BU4 

43  991 

13.8 

-21.8  (105.1) 

-147.8  (103.6) 

1273 

^1.9(176.6) 

+  11.0(171.8) 

BN4 

3095 

14.1 

+4.7  (45.8) 

-229.1  (80.3) 

124 

-68.8  (136.6) 

-11.4  (121.9) 

BT4 

84012 

52.7 

+3.2(107.7) 

-178.9  (111.1) 

1896 

-18.6  (172.2) 

-12.9  (167.0) 

GU6 

76  671 

22.7 

-58.4  (109.7) 

-8.8(91.4) 

1633 

-27.2  (178.3) 

+33.7  (166.2) 

GN6 

4131 

16.3 

-56.3  (51.8) 

-5.4  (87.8) 

150 

+21.8(148.8) 

-11.9(136.5) 

GT6 

93  976 

119.1 

-55.9  (135.4) 

-27.3  (119.2) 

1984 

-5.1  (162.3) 

+0.4(157.3) 

GU4 

77  110 

22.6 

-58.2  (109.9) 

-8.9  (91.4) 

1634 

-26.5  (178.2) 

+33.0(165.5) 

GN4 

4356 

16.2 

-56.6  (51.5) 

+6.2  (86.3) 

154 

+20.1  (150.0) 

-14.0(133.3) 

GT4 

93  978 

119.1 

-55.8  (135.4) 

-27.4  (119.3) 

1983 

-5.2(161.8) 

+0.3  (157.4) 

SU6 

59184 

30.0 

+6.5  (22.3) 

+  14.2  (26.2) 

241 

+5.7  (44.3) 

+25.6  (49.9) 

SN6 

4035 

16.0 

-7.6  (24.2) 

-26.8  (48.9) 

37 

-15.0  (45.7) 

-1.8  (37.7) 

ST6 

74  838 

167.6 

+5.3  (49.4) 

+20.8  (55.5) 

301 

+9.8  (65.9) 

+30.1  (63.1) 

SU4 

59125 

30.0 

+6.4  (21.8) 

+  14.1  (25.8) 

240 

+6.2  (43.3) 

+24.1  (47.6) 

SN4 

4036 

16.0 

-7.5  (23.9) 

-26.5  (48.3) 

37 

-14.7  (45.4) 

-0.8  (37.0) 

ST4 

74  838 

167.6 

+5.4  (49.3) 

+  19.9  (55.2) 

301 

+9.7  (65.9) 

+29.0  (63.3) 

"  Comparative  results  given  by  all  types  of  local  correction  designed.  The  solutions  nomenclature  is  formed  by  one  letter  designating  the  quasar 
position  input  catalog,  followed  by  one  letter  designating  the  stellar  reference  frame,  and  by  the  codes  6  and  4  informing  whether  the  algorithm  is 
of  the  six  or  four-parameter  model.  Nq  is  the  number  of  quasars  for  which  a  given  local  correction  was  possible  and  succeeded.  Ns  is  the  average 
number  of  catalog  stars  contributing  to  the  solution.  Ca  cos  5  and  C6  columns  bring  the  mean  corrections  on  right  ascension  and  declination,  and 
the  standard  deviations  (in  brackets).  The  final  two  columns  show  the  average  optical  to  radio  offset,  preceded  by  the  number  of  sources  entering 
the  average  (Nrad).  All  angular  values  are  in  milli-arcsec.  See  the  text  for  detailed  explanations  on  the  solutions  and  samples. 


plate  where  the  quasar  is  found)  because  of  the  smallness  of  the 
used  field.  The  Bl.O  mean  epoch  is  J1977.12  (+5.49  y).  In  the 
other  hand,  the  Bl.O  positions  are  given  by  an  average  of  plate 
reductions,  in  such  a  way  to  jeopardize  the  hypothesis  of  com¬ 
mon  epochs,  given  the  range  of  magnitudes  and  to  some  extent 
the  two  different  types  of  objects.  For  the  GSC2.3,  in  contrast, 
the  positions  always  refer  to  one  single  highest  quality  plate, 
whenever  possible  the  red  one.  But  then  again,  the  GSC2.3  mean 
epoch  is  J1992.5  (+3.72  y)  and  therefore  there  is  a  balance  be¬ 
tween  reducing  modestly  the  stellar  position  scatter  and  the  er¬ 
rors  that  might  be  introduced  by  faulty  or  biased  proper  motions. 

To  verify  which  way  is  preferable,  the  local  corrections  for 
the  sources  belonging  to  the  enlarged  radio  frame  were  evalu¬ 
ated  both  at  J2000  and  the  quasar  position  epoch.  A  second  trial 
was  made  by  considering  only  the  sources  from  the  enlarged  ra¬ 
dio  frame  that  belonged  to  the  three  input  catalogs  and  could 
be  reduced  by  all  the  stellar  frames.  This  subsample  contains 
240  sources  but  enables  a  direct  comparison  between  all  the  pos¬ 
sible  ways  of  treatment.  The  outcomes  show  no  net  advantage  in 
working  at  the  quasar  epoch.  For  the  Bl.O  catalog,  the  number  of 
sources  for  which  the  local  reduction  converges  drops;  there  is 
a  marginal  gain  on  the  number  of  stars  entering  in  the  reduction 
but  the  number  of  actually  used  stars  remains  the  same;  there  is 
no  decreasing  of  the  standard  deviations  from  the  solution  cor¬ 
rections,  or  from  the  optical  minus  radio  offsets;  and  the  removal 
of  the  equatorial  bias  is  less  efficient  (by  98  mas).  The  external 
check  actually  indicates  a  small  increase  of  the  average  detach¬ 
ment  from  the  representation  of  the  ICRF  origin.  For  the  GSC2.3 
there  is  no  important  loss  in  working  at  the  quasar  epoch,  but  no 
gain  either.  The  number  of  sources  for  which  the  local  reduction 
converges  remains  the  same,  and  so  does  the  number  of  entry 
and  used  stars.  The  standard  deviation  in  the  solution  does  not 
change,  and  for  the  optical  minus  radio  positions  the  offset  scat¬ 
ter  improves  only  marginally  (by  4.5  mas).  The  equatorial  bias 
is  less  efficiently  removed  (by  40  mas).  Exactly  the  same  results 
are  obtained  when  inspecting  the  subsample  of  common  sources. 
In  view  of  these  results,  we  decided  to  work  only  at  J2000  for 


the  local  corrections.  Thus,  all  pairs  of  input  quasar  catalogs  and 
stellar  catalogs  are  treated  identically. 

A  second  point  concerning  the  issue  of  epochs  has  already 
been  touched  in  the  previous  discussion.  It  concerns  whether  the 
local  reductions  made  with  the  UCACN,  that  is  for  the  northern¬ 
most  sources,  and  with  the  2MASS  would  be  significantly  worse 
than  those  made  with  the  UCAC2.  As  we  have  seen,  both  of  them 
have  a  mean  epoch  close  to  J2000,  where  the  local  corrections 
are  performed.  For  the  2MASS  catalog,  which  is  not  originally 
presented  as  an  astrometric  frame,  the  number  of  equations  of 
condition  is  substantially  increased  for  a  same  patch  of  the  skies 
as  defined  for  the  UCAC2.  To  test  the  performance  of  the  differ¬ 
ent  catalogs,  we  used  the  DR5,  which  was  shown  in  Sect.  3  to 
have  intrinsically  smaller  scatter  relative  to  the  ICRF.  As  before, 
the  enlarged  radio  frame  and  its  subset  of  all  common  sources  are 
employed,  and  they  provide  very  similar  responses.  The  UCAC2 
based  reductions  produce  the  best  results  as  expected.  The  right 
ascension  and  declination  standard  deviation  in  the  offsets  to  the 
radio  frame  is  at  47.5  mas.  We  remark  that  this  is  already  a  gain 
relative  to  the  DR5  catalog  offsets  that  are  at  59  mas,  as  produced 
from  Table  1  values.  The  standard  deviation  in  the  solution  it¬ 
self  is  at  33  mas.  The  UCACN  based  reductions  perform  much 
alike,  being  actually  slightly  better  on  declination  and  30  mas 
worse  on  right  ascension.  The  optical  minus  radio  standard  de¬ 
viation  is  14  mas  worse  than  that  of  the  UCAC2.  The  2MASS 
based  reductions  fare  somewhat  worse  than  the  UCACN  ones, 
the  solution  standard  deviation  being  larger  by  18  mas,  and  the 
optical  minus  radio  standard  deviations  being  worse  by  6.5  mas. 
Nonetheless,  we  note  that  the  average  optical  to  radio  position 
offsets  are  much  similar  for  the  UCAC2  and  the  2MASS  (dif¬ 
fering  by  just  5  mas).  Finally,  we  remark  that,  eventhough  not 
bringing  a  direct  improvement  in  accuracy  as  the  UCAC2  based 
reductions,  the  UCACN  and  the  2MASS  local  reductions  enable 
us  to  place  the  DR5  positions  on  the  same  frame  with  the  Bl.O 
and  GSC2.3  reduced  positions,  without  impairing  their  intrinsic 
precision. 

Table  3  summons  up  all  the  applied  local  corrections. 
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Fig.  8.  Equatorial  coordinates  distribution  of  the  local  astrometric  corrections  derived  by  the  six-parameter  model  using  the  UCAC  stellar  reference 
frame  onto  the  USNO  Bl.O  quasar  positions.  The  clear  points  correspond  to  the  quasars  for  which  there  are  ERF  radio  positions. 


As  can  be  seen  in  Table  3,  the  different  quasar  input  catalogs 
(codes  B,  G,  or  S)  require  different  local  corrections,  while  the 
three  stellar  reference  catalogs  (codes  U,  N,  or  T)  deliver  such 
corrections  with  similar  accuracy.  The  local  correction  outputs 
(codes  C)  agree  when  calculated  either  by  the  six-parameter  or 
by  the  four-parameter  models  (codes  6  or  4).  The  local  correc¬ 
tions  are  generally  larger  than  the  remaining  optical  minus  radio 
offsets  (codes  A),  which  supports  their  usefulness.  The  optical 
minus  radio  mean  offset  in  all  cases  is  improved  relatively  to  the 
one  obtained  using  original  quasar  catalogs  positions.  As  a  con¬ 
sequence,  the  differences  between  the  results  from  the  various 
local  reductions  become  similarly  small. 

At  this  point,  the  three  quasar  input  catalogs  and  the  three 
stellar  reference  frames  deliver  concurrent,  independent  sets  of 
positions.  Whenever  the  1st  degree  polynomial  solution  fails 
to  converge,  the  four-parameter  solution  results  can  be  adopted 
equivalently.  There  are  thus  up  to  twelve  possible  local  solu¬ 
tions  calculated  for  each  source.  For  61%  of  the  sources  more 
than  one  solution  was  obtained,  for  844  sources  all  twelve  solu¬ 
tions  were  completed,  and  for  12  sources  no  solution  converged. 
Figures  8  to  1 1  display  the  local  corrections  for  the  most  crit¬ 
ical  cases,  namely  the  Bl.O  and  the  GSC2.3,  as  given  by  the 
UCAC  and  2MASS  frames.  Comparing  with  Figs.  5,  6,  it  is  al¬ 
ready  clear  that  the  local  corrections  bring  the  quasar’s  positions 
onto  the  ICRF.  The  comparison  between  these  hgures  and  the 
initial  plots  (Figs.  5-7)  show  that  the  local  astrometric  solution 
significantly  redress  the  input  positions  from  the  large  catalogs. 


We  note  that  the  outcomes  from  the  UCAC  and  2MASS  stellar 
reference  frames  have  much  in  common.  The  six-parameter  and 
the  four-parameter  models  also  provide  coincident  corrections, 
and  thus,  only  those  from  the  six-parameter  model  are  shown  in 
the  hgures.  We  also  note  that  there  is  no  discontinuity  between 
the  zones  where  either  the  UCAC2  or  the  UCACN  stellar  frames 
were  used. 


4.  The  global  orientation  towards  the  ICRS 

In  this  section,  points  7  and  8  presented  in  Sect.  2.11  (data  How) 
are  derived. 

Once  the  local  correction  is  obtained,  the  next  move  is  to 
tie  the  preliminary  frame  just  derived  to  the  ICRF.  As  discussed 
in  Sect.  1,  the  ICRF  alone  is  too  sparse,  and  in  particular  the 
number  of  its  bright  optical  counterparts,  to  suffice  for  a  detailed 
tying.  Since  the  long  base  radio  interferometric  observations  al¬ 
ways  combine  data  for  several  sources  in  different  sessions,  a 
radio  position  net  that  is  considerably  larger  than  the  ICRF  can 
be  formed  by  adding  the  VLBA  calibrators  and  the  VLA  cali¬ 
brators.  In  the  case  of  the  VLA,  only  sources  of  the  highest  as¬ 
trometric  precision  are  included.  The  accurate,  ICRF  represent¬ 
ing,  radio  interferometric  positions  of  the  enlarged  radio  sample 
sources  (ERF)  are  those  compared  with  their  optical  counterpart 
from  the  stellar  based  preliminary  frame  derived  from  the  lo¬ 
cal  corrections.  From  the  ICRF,  718  sources  are  collected,  while 
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Fig.  9.  Equatorial  coordinates  distribution  of  the  local  astrometric  corrections  derived  by  the  six-parameter  model  using  the  2MASS  stellar  refer¬ 
ence  frame  onto  the  USNO  Bl.O  quasar  positions.  The  clear  points  correspond  to  the  quasars  for  which  there  are  ERF  radio  positions. 


from  the  VSC6  and  the  VLAC,  2684  sources  and  59  sources  are 
collected  respectively. 

The  enlarged  radio  frame  produced  in  this  way  thus  contains 
3461  sources.  They  are  adequately  represented  in  all  of  the  three 
contributing  catalogs  as  described  in  Table  3  (Col.  Nrad).  The 
sky  distribution  was  shown  in  Fig.  4.  From  the  ERF  sources, 
2263  have  an  optical  counterpart  in  at  least  one  of  the  input  cat¬ 
alogs.  The  average  minimum  distance  between  the  ERF  sources 
is  187  arcmin,  with  a  mode  of  70  arcmin. 

The  ERF  radio  positions  are  then  compared  with  those  of 
their  optical  counterparts  from  the  stellar  based  preliminary 
frame  derived  by  the  local  corrections.  Initially,  the  compari¬ 
son  entails  the  derivation  of  global  orientation  corrections  about 
the  equatorial  coordinates  standard  triad,  plus  a  correction  for 
a  bias  relatively  to  the  equator.  Following  (Arias  et  al.  1988), 
the  canonical  equations  relating  the  right  ascension  and  declina¬ 
tion  offsets  to  the  rotations  (A1-A3)  about  the  equatorial  axes 
(namely,  X  on  the  equator  pointing  to  the  conventional  origin, 
Y  on  the  equator  perpendicular  to  X,  and  Z  pointing  to  the  con¬ 
ventional  pole)  are  given  by 

Aa  cos  (5  =  Al  sin  d  cos  a  +  A2  sin  6  sin  a  -  A3  cos  6  (5) 

Ad  =  -Al  sin  a  +  A2  cos  a  -i-  A4.  (6) 

Equations  (5)  and  (6)  are  solved  independently  to  assess  the  ro¬ 
bustness  of  the  corrections.  In  this  case,  the  direction  cosines  A I 
and  A2  agree  when  calculated  either  from  Aa  cos  d  or  Ad  within 


1  cr  in  most  cases.  Eigure  12  shows  the  good  agreement  on  Al 
and  A2  as  calculated  from  the  right  ascension  residuals,  from 
the  declination  residuals,  or  from  the  combined  residuals.  The 
equatorial  bias  (which  is  discussed  in  the  next  section)  can  also 
be  independently  derived  as  simply  Ad  =  A4.  Also  for  A4,  either 
the  straight  average  or  the  adopted  combined  solution  agree  well. 

Eor  the  sake  of  gathering  the  largest  number  of  equations 
of  condition,  the  hnal  values  are  derived  by  a  combined  matrix 
including  the  three  rotation  corrections,  plus  the  representation 
of  the  equatorial  bias. 

4.1.  The  equatorial  bias 

Former  investigations  (Assahn  et  al.  2001)  pointed  out  marked 
differences  in  the  equatorial  bias  concerning  the  north  and  south 
portions  of  the  USNO  A2.0  catalog.  The  point  of  discontinuity 
is  around  -20°  of  declination,  which  is  the  boundary  of  the  north 
and  south  plate  surveys  that  contributed  to  that  catalog.  Although 
such  marked  discontinuity  should  not  befall  neither  the  Bl.O  nor 
the  GSC23  catalogs,  they  mostly  share  the  same  constituent  sur¬ 
veys  with  the  USNO  A2.0. 

To  analyze  whether  a  discontinuity  remains  in  the  data, 
Eq.  (6)  was  solved  for  the  catalogs  Bl.O  and  GSC23,  as  locally 
corrected  by  the  UCAC2  and  the  2MASS  frames,  by  adopting 
a  single  A4  unknown,  two  A4  unknowns  dividing  the  Ad  off¬ 
sets  at  -20°  of  declination,  and  eighteen  A4  unknowns  to  fit  the 
Ad  offsets  within  declination  strips  10°  wide.  The  plots  in  Eig.  13 
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Fig.  10.  Equatorial  coordinates  distribution  of  the  local  astrometric  corrections  derived  hy  the  six-parameter  model  using  the  UCAC  stellar  refer¬ 
ence  frame  onto  the  GSC2.3  quasar  positions.  The  clear  points  correspond  to  the  quasars  for  which  there  are  ERF  radio  positions. 


show  features  that  are  more  complex  than  could  he  described  by 
one  A4  standard  offset.  There  is  a  coarse  positive  plateau  south¬ 
wards  of  declination  -20°,  and  values  averaging  near  zero  north¬ 
wards.  The  standard  deviation  in  the  residuals  from  the  18  strips 
solutions  is  indeed  no  smaller  than  that  of  the  -20°  boundary, 
2  A4  terms  solution.  This  is  indicative  of  an  interplay  of  fac¬ 
tors  that  can  be  more  adequately  handled  in  the  next  step  where 
harmonic  terms  are  introduced.  On  the  other  hand,  the  solution 
contemplating  a  single  A4  term  fares  the  worst  by  a  few  mas, 
and  is  certain  to  provide  the  least  appropriate  description. 

The  A4  declination  dependency  that  is  verified,  prompted 
us  to  assess  whether  a  similar  right  ascension  dependency  is 
present.  In  this  case,  no  such  dependency  exists. 

4.2.  Rotation  angles  and  equatorial  bias 

The  final  values  of  the  Al,  A2,  and  A3  direction  cosines  for  the 
global  orientation  towards  the  ICRF,  as  well  as  the  equatorial 
bias,  were  found  by  adjusting  the  optical  minus  radio  residual 
between  the  positions  evaluated  from  the  local  astrometric  solu¬ 
tion  and  those  given  by  the  enlarged  radio  frame.  The  values  are 
presented  in  Table  4. 

For  the  Bl.O  and  GSC2.3  northernmost  sources,  for  which 
the  local  correction  was  made  using  the  UCACN,  separate  val¬ 
ues  were  derived  for  the  rotation  angle  and  the  (northern  only) 
equatorial  bias.  But,  for  the  DR5  northernmost  sources,  al¬ 
though  the  number  of  sources  handled  with  UCACN  is  small,  the 


precisions  are  high  and  so  a  complete  solution  including  only 
one  equatorial  bias  term  was  sufficient.  We  remark  that,  for  the 
DR5,  the  Al  and  A2  values,  although  significant  and  represent¬ 
ing  true  corrections  to  be  applied  to  the  quasar  equatorial  coor¬ 
dinates,  must  be  interpreted  with  caution  in  terms  of  what  they 
would  represent  for  the  DR5  frame  orientation,  because  of  the 
concentration  of  the  DR5  sources  in  the  direction  of  the  right 
ascension  of  12h. 


5.  The  removal  of  zonal  bias 

In  this  section,  points  9  and  10  of  Sect.  2.11  (data  flow)  are 
derived. 

5.1.  Harmonic  functions 

At  this  point,  the  optical  frame  was  homogenized  by  the  local 
corrections,  and  was  globally  aligned  with  the  ICRF.  The  fi¬ 
nal  point  is  to  sweep  its  axes  to  check  for  any  remaining  zonal 
modulations.  This  was  done  for  the  right  ascension,  declination, 
and  R  magnitude  dimensions.  In  a  sense,  this  is  equivalent  to 
straighten  up  the  axes  from  waves  of  corrugation  either  trapped 
in  the  stellar  catalogs  themselves,  or  that  were  present  in  the  ini¬ 
tial  deep  catalogs  and  could  not  be  completely  removed  by  the 
local  corrections. 

The  problem,  the  quantity  of  assessing  points,  the  sky  cov¬ 
erage,  the  dimensions  of  dependencies,  and  the  relationship 
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Fig.  11.  Equatorial  coordinates  distribution  of  the  local  astrometric  corrections  derived  by  the  six-parameter  model  using  the  2MASS  stellar 
reference  frame  onto  the  GSC2.3  quasar  positions.  The  clear  points  correspond  to  the  quasars  for  which  there  are  ERF  radio  positions. 


between  the  less  precise,  corrugated  frame  and  the  more  precise, 
correcting  one,  are  all  the  same  in  the  case  considered  by  Schwan 
(1988)  when  deriving  zonal  corrections  to  the  FK5  construction. 
Thus,  the  same  treatment  by  spherical  harmonics  (Brosche  1966) 
is  used  here.  The  right  ascension  and  declination  corrections  are 
expressed  by  the  linear  combination  of  the  terms  of  a  series,  in 
which  the  cpnml-ih  element  is  given  by 

A  c  =  Ccpnml  Hp(R')  Ln(sm  6)  Fml(a)  (7) 

where  Ac  is  one  or  the  other  equatorial  coordinates  offsets,  since 
independent  spherical  harmonic  series  are  derived  for  each  of 
them;  Ccpnml  is  the  cprml-th  coefficient  to  be  derived,  for  the 
c  equatorial  coordinate  series;  Hp(R')  is  the  Hermite  polyno¬ 
mial  of  degree  p,  which  is  a  function  of  the  R  magnitude  nor¬ 
malized  by  the  mean  and  standard  deviation  for  the  quasars; 
Ln(sin  6)  is  the  Legendre  polynomial  of  degree  n,  function  of 
sin  (5;  and  Fml(a)  is  the  Fourier  term  of  order  mona  and  sign  1. 

The  offsets  are  produced  by  comparing  the  equatorial  coor¬ 
dinates  from  the  previous  step  (global  rotation)  with  the  VLBI 
positions  from  the  ERF.  That  is,  2.263  pairs  are  available.  This 
is  again  similar  to  the  number  that  were  available  for  the  Schwan 
(1988)  application.  Therefore,  terms  up  to  very  higher  order  can 
be  derived.  However,  the  analysis  from  the  earlier  steps  advises 
us  to  be  cautious  in  defining  the  significant  order  up  which  to 
halt.  The  distribution  in  right  ascension  of  the  offsets  from  the 
original  deep  catalogs  suggests  a  In  term  and  smaller  ripples 
(cf.  Figs.  5-7),  while  the  analysis  of  global  rotation  coefficients 


Solution 

Fig.  12.  Coincidence  of  the  Al  and  A2  direction  cosines  as  obtained 
either  from  the  Aa  cos  S  residuals,  from  the  AS  residuals,  or  from  both. 
The  value  adopted  was  from  the  solution  including  both  the  right  ascen¬ 
sion  and  declination  residuals. 


exhibits  no  a  dependence.  In  contrast,  the  analysis  of  the  equa¬ 
torial  bias  displays  a  6  dependence,  coherent  in  bins  of  10° 
(Fig.  13).  Magnitude  dependences  are  believed  to  be  possible 
within  the  large  quasar  catalogs.  They  could  to  some  extent  be 
introduced  by  the  local  solution  itself,  because  of  the  magni¬ 
tude  difference  between  stars  and  quasars.  The  analysis  made 
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Fig.  13.  Declination  dependency  of  the  equatorial  bias  (A4).  The  plotted  averages  correspond  to  hins  10°  wide  in  declination.  The  shaded  codes 
indicate  the  numher  of  ERF  quasars  within  the  hins. 


in  Sect.  3.3  points  out  to  such  dependences  but  of  small  ampli¬ 
tude.  All  accounted  for,  and  to  have  representative  quantities  of 
ERF  objects  with  which  to  compute  the  coefficients  of  each  term, 
we  adopted  to  extend  to  9th  order  the  Lagrange  (sin  6)  polyno¬ 
mials  and  the  Fourier  (or)  functions.  The  Hermite  {R  magnitude) 
polynomials  were  extended  to  2nd  order,  but  just  to  5th  order 
of  the  other  two  terms.  Given  the  poor  sky  distribution,  declina¬ 
tion  terms  were  not  used  for  catalog  combinations  including  the 
UCACN  frame,  nor  right  ascension  terms  for  catalog  combina¬ 
tions  including  the  DR5  quasar  list.  All  the  above  choices  were 
made  to  retain  only  truly  significant  terms  derived  from  the  ERF 
comparison,  such  that  they  could  be  confidently  applied  to  the 
full  set  of  quasars. 

Since  by  definition  the  spherical  harmonics  are  orthogonally 
independent,  their  statistical  significance  can  be  tested  individu¬ 
ally.  Each  term  was  initially  tested  separately.  Only  those  statis¬ 
tically  significant  at  3cr  with  coefficients  not  smaller  than  1  mas 
were  retained.  Once  more,  this  taxes  heavier  on  the  high  order 
polynomials.  Next,  starting  from  the  lower  orders  and  up,  the 
corrections  brought  by  each  term  are  removed  from  the  equa¬ 
torial  offsets,  and  the  following  term  is  individually  tested.  The 
same  criterion  of  3cr  significance  and  coefficient  not  smaller  than 
1  mas  is  applied  to  retain  the  terms.  At  this  round,  the  high¬ 
est  order  Hermite  magnitude  polynomials,  when  convolved  with 
the  higher  Legendre  (declination)  and  Fourier  (right  ascension) 
orders,  are  retained  only  if  the  coefficient  is  not  smaller  than 
5  mas.  The  terms  surviving  these  two  rounds  are  added  by  a 
linear  combination  to  adjust  the  right  ascension  and  declination 
offsets,  separately,  by  least  squares.  Table  5  reports  the  number 
of  terms  surviving  the  rounds.  We  note  the  important  decreasing 


rate  from  the  first  to  the  second  round,  which  justifies  the  rounds 
scheme.  On  the  other  hand,  there  is  no  dropping  from  the  second 
round  (one  by  one  regimen)  to  the  third  round  (least  squares  fit 
of  all  validated  terms  together),  what  shows  that  they  are  truly 
statistically  significant. 

Table  6  reports  the  coefficients  derived  during  the  final  ad¬ 
justment  for  each  family  of  quasar  catalog  and  stellar  frame.  As 
expected,  the  coefficients  do  not  vary  from  the  second  round  to 
the  final  adjustment,  indicating  that  no  secondary  harmonic  sur¬ 
vived.  There  is  only  one  significant  term  including  the  DR5.  It 
relates  to  the  2MASS  local  correction  and  depends  solely  on  dec¬ 
lination  (3rd  order).  This  reflects  the  high  quality  of  the  DR5  po¬ 
sitions  and  indicates  that  the  local  correction  did  not  induce  any 
derail  on  the  input  quasar  catalogs.  For  the  UCACN  there  is  also 
only  one  significant  term.  It  relates  to  the  GSC2.3  quasar  catalog 
and  is  of  high  right  ascension  order  (8th).  The  lack  of  UCACN 
harmonic  terms  reflects  the  confinement  of  the  sky  region  where 
it  was  applied.  The  positions  that  originated  in  the  Bl.O  required 
42  harmonic  terms  (3 1  for  the  right  ascension  offsets),  and  those 
from  the  GSC23  32  harmonic  terms  (20  for  the  right  ascension 
offsets).  Considering  the  stellar  reference  frames,  the  UCAC2 
required  32  harmonic  terms  (25  for  the  right  ascension  offsets), 
while  the  2MASS  required  42  harmonic  terms  (28  for  the  right 
ascension  offsets).  In  all  cases,  the  right  ascension  offsets  re¬ 
quired  more  harmonic  terms.  In  contrast,  the  Ln(sind)  domi¬ 
nated  terms  in  general  correspond  to  larger  coefficients,  which 
worked  to  correct  the  deviations  detected  at  the  A4  equatorial 
bias  analysis.  The  magnitude  terms  appear  only  at  the  1st  or¬ 
der,  in  18  Bl.O  related  terms  and  10  GSC23  related  terms.  Thus, 
magnitude  terms  are  presented  as  the  catalogs  authors  cautioned. 
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Table  4.  Direction  cosines  relative  to  the  ERF  after  the  application  of  the  local  astrometric  corrections". 


Solution 

N 

Al 

A2 

A3 

A4 

A5 

BU6 

2616 

-33.7  ±  6.3 

-8.0  ±6.1 

50.8  ±  5.3 

40.5  ±  10.7 

2.6  ±5.3 

BN6 

244 

-5.5  ±  12.9 

31.8  ±  14.1 

123.0  ±25.2 

1.5  ±  12.8 

BT6 

3764 

-23.8  ±5.1 

1.5  ±4.7 

23.7  ±  4.6 

27.5  ±  9.7 

-21.0  ±4.3 

BU4 

2522 

-37.5  ±  6.3 

-4.9  ±  6.2 

48.6  ±  5.4 

42.1  ±  10.7 

0.0  ±  5.4 

BN4 

240 

-0.4  ±  12.9 

34.7  ±  14.1 

118.9  ±  25.4 

-0.3  ±  12.8 

BT4 

3762 

-24.6  ±5.1 

1.4  ±4.7 

24.3  ±  4.6 

24.7  ±  9.6 

-21.7  ±4.3 

GU6 

3238 

-7.0  ±  5.8 

-3.1  ±5.3 

30.4  ±  4.7 

100.8  ±  9.3 

15.2  ±4.8 

GN6 

300 

-7.2  ±  13.2 

16.3  ±  14.2 

-43.4  ±25.1 

-4.7  ±  13.1 

GT6 

3940 

^.0±4.7 

-4.9  ±  4.4 

8.7  ±4.3 

30.2  ±  8.7 

-6.1  ±4.0 

GU4 

3240 

-7.7  ±  5.8 

-2.5  ±  5.3 

29.4  ±  4.7 

99.1  ±9.3 

14.7  ±4.8 

GN4 

308 

-7.1  ±  13.1 

15.7  ±  14.0 

-40.4  ±  24.7 

-6.8  ±  12.9 

GT4 

3938 

-3.7  ±  4.7 

-6.0  ±  4.4 

8.7  ±4.3 

30.1  ±8.7 

-6.3  ±  4.0 

SU6 

480 

-5.3  ±  4.7 

-14.6  ±  4.8 

^.8  ±  3.7 

16.6  ±4.1 

SN6 

74 

-5.0  ±  11.6 

14.8  ±  13.9 

32.0  ±  17.8 

9.8  ±  13.1 

ST6 

600 

-13.9  ±5.5 

-14.6  ±  6.0 

-5.1  ±4.9 

20.1  ±5.3 

SU4 

478 

-5.5  ±  4.5 

-12.8  ±4.6 

-5.2  ±  3.5 

16.3  ±4.0 

SN4 

74 

-6.2  ±  11.4 

13.1  ±  13.8 

33.2  ±  17.6 

9.5  ±  13.0 

ST4 

600 

-14.3  ±  5.5 

-14.2  ±  6.0 

-4.7  ±  4.9 

19.2  ±5.3 

"  Direction  cosines  Al,  A2,  and  A3  from  the  comparison  of  the  equatorial  positions  as  found  for  each  family  of  solutions  and  as  given  by  the  radio 
interferometry  for  the  ERF  objects.  The  equatorial  bias  A4  (north  of  declination  -20°)  and  A5  (south  of  declination  -20°),  when  applicable,  are 
also  presented.  The  solutions  nomenclature  is  formed  by  one  letter  designating  the  quasar  position  input  catalog,  followed  by  one  letter  designating 
the  stellar  reference  frame,  and  by  the  codes  6  and  4  informing  whether  the  algorithm  is  of  the  six  or  four-parameter  models.  N  is  the  number  of 
equations  used  to  derive  the  direction  cosines  and  equatorial  biases.  All  angular  values  are  in  milli-arcsec. 


but  they  are  less  important  than  the  zonal  terms.  Figure  14  shows 
the  histograms  of  the  harmonic  terms  coefficients  and  associated 
errors.  The  errors  resemble  closely  a  Poisson  distribution.  The 
coefficients  peak  around  -6  mas  for  the  right  ascension  offsets, 
and  around  -t-6  mas  for  the  declination  offsets. 

5.2.  Local  inhomogeneities 

Figures  15  and  16  (open  symbols)  show  that  after  the  correction 
by  harmonic  series,  the  optical  minus  radio  residuals  are  equal  to 
zero  within  the  statistical  significance  of  3  x  standard  error,  that 
is  without  being  attached  to  the  scatter  of  the  optical  position 
determinations.  This  is  shown  for  the  pairs  combining  the  B 1 .0 
(Fig.  15)  and  GSC  (Fig.  16)  quasar  input  lists,  and  the  UCAC2 
(squares)  and  2MASS  (circles)  local  stellar  frames,  which  are 
the  largest  and  have  the  most  homogenous  sky  coverage  of  the 
combinations  used  here. 

Notwithstanding  this,  some  clumping  is  seen,  without  a  clear 
pattern  over  the  sky,  that  can  be  dealt  with  further.  The  ERF  sky 
density  is  such  that  the  mean  distance  from  any  given  quasar  to 
the  closest  ERF  VLBI  position  is  1.7°.  And,  on  average,  there  are 
10  ERF  VLBI  positions  within  a  radius  of  4.4°  around  any  given 
quasar,  the  largest  distance  being  on  average  6.3°.  Under  such 
conditions  for  each  of  the  combinations  of  quasar  list  and  stellar 
reference  frame,  hence  for  every  quasar  a  correction  for  the  lo¬ 
cal  inhomogeneities  towards  the  VLBI  position  can  be  obtained 
from  the  average  of  10  ERF  quasars.  Within  the  ERF  quasars 
radius,  objects  are  removed  one  by  one  if  both  their  right  as¬ 
cension  and  declination  offsets  are  larger  than  2cr.  On  average, 
1.3  ERF  quasars  are  removed,  what  shows  the  robustness  of  the 
local  inhomegeneities  corrections.  This  procedure  is  similar  to 
that  of  overlapping  circles  (Taff  et  al.  1990a,b),  but  no  over¬ 
lap  is  actually  required  since  the  autocorrelation  between  adja¬ 
cent  corrections  fades  into  a  characteristic  value  at  a  distance 
of  20°.  Figure  17  illustrates  the  weakness  of  the  autocorrelations 
by  taking  each  quasar  in  turn  as  a  pole  source,  and  plotting  the 


Coefficients  (mas)  Errors  (mas) 

Fig.  14.  Distribution  of  the  coefficients  and  associated  errors  of  the  har¬ 
monic  functions  fitting  the  Aa  cos  6  and  the  AS  residuals.  All  significant 
terms  (first  round  -  see  text)  are  counted  here. 


standard  deviation  in  the  inhomogeinity  correction  within  grow¬ 
ing  distance  rings.  For  nearby  sources,  i.e.,  very  small  rings,  the 
inhomogeneities  are  much  the  same,  and  the  standard  deviation 
is  accordingly  small.  Rapidly,  however,  before  reaching  the  dis¬ 
tance  of  20°,  the  standard  deviation  already  reaches  the  limit  of 
scatter  characteristic  of  each  family  of  input  catalog  and  stellar 
reference  frame  solutions. 

By  determining  in  this  way  the  correction  for  the  local  inho¬ 
mogeneities,  it  is  seen  in  Figs.  15  and  16  (filled  symbols)  that  the 
clumps  are  clearly  minimized.  After  applying  the  corrections  by 
harmonic  functions,  the  optical  minus  radio  offsets  are  typically 
null  within  1.3cr,  what  is  reduced  further  to  within  0.7cr  after 
applying  the  complementary  local  inhomogeneities  corrections. 
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Table  5.  Number  of  significant  harmonic  terms  in  the  three  rounds  of  testing". 


Solution 

RA 

Dec 

1st  round 

2nd  round 

3rd  round 

1st  round 

2nd  round 

3rd  round 

BU6 

21 

6 

6 

2 

1 

1 

BN6 

0 

0 

0 

0 

0 

0 

BT6 

17 

9 

9 

17 

5 

5 

BU4 

23 

7 

7 

0 

0 

0 

BN4 

0 

0 

0 

0 

0 

0 

BT4 

17 

9 

9 

17 

5 

5 

GU6 

15 

6 

6 

8 

3 

3 

GN6 

0 

0 

0 

2 

1 

1 

GT6 

9 

4 

4 

5 

2 

2 

GU4 

15 

6 

6 

8 

2 

2 

GN4 

0 

0 

0 

2 

1 

1 

GT4 

9 

4 

4 

5 

2 

2 

SU6 

0 

0 

0 

0 

0 

0 

SN6 

0 

0 

0 

0 

0 

0 

ST6 

1 

1 

1 

0 

0 

0 

SU4 

0 

0 

0 

0 

0 

0 

SN4 

0 

0 

0 

0 

0 

0 

ST4 

1 

1 

1 

0 

0 

0 

"  Harmonic  terms  surviving  in  the  three  rounds  of  testing  for  the  most  significant  ones.  The  three  rounds  are  explained  in  the  text,  and  the  third 
round  corresponds  to  the  adopted  harmonic  fit.  All  procedure  was  done  independently  for  the  right  ascension  and  declination  residuals  between 
the  solutions  and  the  ERF  radio  interferometry  positions.  The  solutions  nomenclature  is  formed  by  one  letter  designating  the  quasar  position  input 
catalog,  followed  by  one  letter  designating  the  stellar  reference  frame,  and  by  the  codes  6  and  4  informing  whether  the  algorithm  is  of  the  six  or 
four-parameter  model. 
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Fig.  15.  Averaged  Aor  cos  S  and  AS  optical  (from  the  USNO  B 1 .0)  minus  radio  (from  the  ERF)  residuals  in  units  of  their  standard  deviations.  Open 
symbols  refer  to  the  values  corrected  by  the  harmonic  functions.  Filled  symbols  refer  to  the  values  further  redressed  by  the  local  inhomogeneities 
corrections.  In  both  cases  the  averages  are  non  significant  below  3cr.  The  actual  significance  averages  are  1.2cr  after  the  harmonics  treatment  and 
0.8(T  after  the  complementary  local  inhomogeneities  corrections. 
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Fig.  16.  Averaged  Aacos^  and  AS  optical  (from  the  GSC2.3)  minus  radio  (from  the  ERF)  residuals  in  units  of  their  standard  deviations.  Open 
symbols  refer  to  the  values  corrected  by  the  harmonic  functions.  Filled  symbols  refer  to  the  values  further  redressed  by  the  local  inhomogeneities 
corrections.  In  both  cases  the  averages  are  non  significant  below  3cr.  The  actual  significance  averages  are  1.3cr  after  the  harmonics  treatment 
and  0.7cr  after  the  complementary  local  inhomogeneities  corrections. 


6.  The  LQRF  catalog 

The  LQRF  catalog  right  ascension  and  declination  equatorial  co¬ 
ordinates  are  calculated  from  the  values  obtained  in  the  final  step 
of  the  data  treatment  (removal  of  zonal  bias)  by  a  weighted  aver¬ 
age.  In  the  average,  for  each  quasar,  the  solutions  from  the  three 
input  lists  (Bl.O,  GSC23,  and  DR5)  and  the  three  stellar  refer¬ 
ence  frames  (UCAC2,  UCACN,  and  2MASS)  are  combined.  The 
weights  are  the  inverse  of  the  square  root  sum  of  the  internal  and 
external  errors  at  each  of  the  correction  steps  plus  the  formal  er¬ 
rors  of  each  entry  quasar  input  list.  To  not  carry  too  high  or  too 
low  a  weight  for  any  given  combination,  however,  for  each  of 
them  separately  the  errors  have  been  assigned  by  quartiles.  The 
solutions  derived  with  different  weighting  schemes  or  even  with¬ 
out  any  weighting  agree  to  within  10  mas.  This  indicates  that  the 
goal  of  homogenizing  the  different  inputs  was  accomplished. 

Since  the  final  positions  originate  in  the  average  of  multiple 
combinations,  attaching  the  final  error  of  them  to  the  squared 
sum  of  the  contributing  errors  would  favor  the  sources  with 
fewer  solutions  obtained.  Likewise,  if  attaching  the  error  to  the 
dispersion  in  the  multiple  contributing  solutions.  It  is  more  real¬ 
istic  to  derive  the  final  error  in  the  equatorial  coordinates  of  each 
source  from  the  external  comparison  of  the  LQRF  position  with 
the  radio  position.  Since  in  the  construction  of  the  LQRF  the 
three  correction  steps  focused  on  the  achievement  of  an  homo¬ 
geneous  frame,  this  is  done  by  taking  the  average  of  the  LQRF 


minus  radio  positions  for  the  10  closest  ERF  sources  around 
each  quasar.  As  seen  before,  this  implies  a  typical  radius  of  4.4° 
around  any  LQRF  quasar.  Figure  18  presents  the  histograms  of 
the  errors  in  the  coordinates.  The  distributions  can  be  identified 
as  Poisson  disfribufions  peaking  af  139  mas  in  right  ascension 
and  at  130  mas  in  declination. 

The  LQRF  contains  100165  sources.  On  average,  every 
source  has  a  neighbor  within  10  arcmin.  Figure  19  shows  the 
sky  distribution  counts  within  squares  of  10°.  Empty  boxes  are 
found  only  in  the  Galactic  plane  towards  the  southern  hemi¬ 
sphere.  Eigure  20  presents  a  vector  map  of  the  offsets  versus  the 
VLBI  positions,  in  bins  of  side  30°  and  at  least  5  quasars.  The 
mean  offset  is  32.7  mas  and  their  distribution  appears  as  ran¬ 
dom,  with  the  largest  offsets  corresponding  to  bins  of  only  5  or 
6  quasars. 

Eigures  21-24  show  the  optical  minus  radio  coordinates 
offsets.  A  comparison  with  the  equivalent  Eigs.  5-7  illustrates 
clearly  the  evenness  of  the  ICRF  representation  produced  by  the 
LQRF  positions. 

The  LQRF  catalog  is  available  by  means  of  CDS  access  (see 
URL  in  the  first  page  footnote).  An  extract  of  its  first  page  ap¬ 
pears  in  Fig.  25.  Its  columns  contain 

1.  the  J2000  right  ascension  (h,  m,  s); 

2.  the  J2000  declination  (d, ',  "); 

3.  the  external  right  ascension  error  (x  cos  6)  (mas); 
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Table  6.  Values  of  the  coefficients  of  the  adjusted  harmonic  functions". 


Solution 

Offset 

Hp(R') 

Ln(sin  d) 

Fml{a) 

Coefficient 

BU6 

Aa  cos  6 

0 

0 

1 

-1 

-h18.7±5.5 

BU6 

Aa  cos  6 

0 

0 

1 

+  1 

+41.7  ±  5.2 

BU6 

Aa  cos  6 

0 

4 

0 

0 

-21.0  +  6.0 

BU6 

Aa  cos  6 

1 

0 

0 

0 

-5.5  ±2.7 

BU6 

Aa  cos  6 

1 

3 

0 

0 

+13.1  ±  2.9 

BU6 

Aa  cos  6 

1 

4 

1 

-1 

+  8.9  ±  2.8 

BU6 

A6 

0 

8 

3 

-Hi 

-20.2  ±  6.0 

BN6 

Aa  cos  6 

_ 

_ 

_ 

_ 

BN6 

AS 

- 

- 

- 

- 

BT6 

Aa  cos  6 

0 

0 

1 

-1 

+16.6  +  4.7 

BT6 

Aa  cos  6 

0 

0 

1 

-Hi 

+38.0  ±4.2 

BT6 

Aa  cos  6 

0 

3 

1 

-Hi 

+  12.3  +  4.4 

BT6 

Aa  cos  6 

0 

4 

0 

0 

-15.1+4.6 

BT6 

Aa  cos  6 

0 

4 

4 

-1 

-18.2  +  4.5 

BT6 

Aa  cos  6 

1 

0 

0 

0 

-8.7 +  2.1 

BT6 

Aa  cos  6 

1 

3 

0 

0 

+  7.5 +  2.1 

BT6 

Aa  cos  6 

1 

5 

0 

0 

-7.0 +  2.1 

BT6 

Aa  cos  6 

1 

5 

1 

-Hi 

+  7.2  +  2.2 

BT6 

AS 

0 

6 

1 

-1 

+  14.9  +  4.8 

BT6 

AS 

0 

8 

3 

-Hi 

-17.3  +  4.5 

BT6 

AS 

0 

9 

6 

-Hi 

-14.2  +  4.6 

BT6 

AS 

1 

3 

1 

-Hi 

+  6.9  +  2.0 

BT6 

AS 

1 

3 

3 

-1 

+  7.2 +  2.1 

GU6 

Aa  cos  6 

0 

0 

1 

+ 

+  18.0  +  4.7 

GU6 

Aa  cos  6 

0 

1 

0 

0 

+64.2  +  9.0 

GU6 

Aa  cos  8 

0 

3 

0 

0 

-36.3  +  8.5 

GU6 

Aa  cos  8 

0 

5 

0 

0 

+22.0  +  6.1 

GU6 

Aa  cos  6 

0 

9 

0 

0 

+  19.9  +  5.4 

GU6 

Aa  cos  6 

1 

0 

1 

+ 

-7.1  +2.4 

GU6 

AS 

0 

1 

0 

0 

-23.0  +  6.1 

GU6 

AS 

0 

2 

4 

-Hi 

-18.6  +  5.7 

GU6 

A6 

1 

0 

0 

0 

+  7.3  ±  2.3 

GN6 

Aa  cos  8 

_ 

_ 

_ 

_ 

GN6 

AS 

0 

0 

9 

-1 

-33.9+  1,1 

GT6 

Aa  cos  6 

0 

0 

1 

-Hi 

+22.4  +  4.0 

GT6 

Aa  cos  8 

0 

1 

0 

0 

+26.0  +  4.6 

GT6 

Aa  cos  8 

0 

3 

0 

0 

-21.0  +  4.4 

GT6 

Aa  cos  6 

1 

1 

1 

-1 

-7.9  ±2.2 

GT6 

AS 

1 

0 

0 

0 

+  6.8  ±  2.1 

GT6 

AS 

1 

2 

4 

-Hi 

+  7.0  +  2.1 

DU6 

Aa  cos  S 

_ 

_ 

_ 

_ 

DU6 

AS 

- 

- 

- 

- 

DN6 

Aa  cos  S 

_ 

_ 

_ 

_ 

DN6 

AS 

- 

- 

- 

- 

DT6 

Aa  cos  S 

0 

4 

0 

0 

-14.8  +  3.8 

DT6 

AS 

— 

— 

— 

— 

"  As  the  results  from  the  six-parameter  and  four-constants  local  as- 
trometreic  solutions  are  quite  alike,  only  the  results  from  the  six- 
parameter  solutions  are  displayed.  The  code  for  the  solutions  are  as  in 
the  previous  tables  and  the  coefficients  are  in  milli-arcsec. 


Fig.  17.  Fast  fading  of  the  autocorrelation  of  the  local  inhomogeneities 
correction.  In  the  figure  is  plotted  the  standard  deviation  of  the  correc¬ 
tions  in  samples  taken  in  rings  of  increasing  radius.  The  upper  lines 
correspond  to  the  USNO  Bl.O  and  GSC2.3  families  of  solutions,  while 
the  lower  lines  correspond  the  SDSS  DR5  families  of  solutions. 


LQRF  external  R.A.  errors  (mas) 


LQRF  external  DEC.  errors  (mas) 


Fig.  18.  Distribution  of  the  external  errors  assigned  to  the  right  ascen¬ 
sion  and  declination  J2000  positions  of  the  LQRF  sources.  The  distri¬ 
butions  can  be  reckoned  to  Poisson  distributions  peaking  at  139  mas  on 
right  ascension  and  at  130  mas  on  declination. 

8.  the  right  ascension  offsets  to  the  radio  interferometry  posi¬ 
tion  (blank  if  there  is  no  radio  position)  (x  cos  6)  (mas); 

9.  the  declination  offsets  to  the  radio  interferometry  position 
(blank  if  there  is  no  radio  position)  (mas); 

10.  the  J2000  radio  interferometry  right  ascension  (h,  m,  s) 
(blank  if  there  is  no  radio  position); 

1 1 .  the  J2000  radio  interferometry  declination  (h,  m,  s)  (blank  if 
there  is  no  radio  position). 


4.  the  external  declination  error  (mas); 

5.  the  number  of  radio  interferometry  neighbor  positions  used 
to  evaluate  the  external  errors; 

6.  the  R  magnitude,  taken  from  the  LQAC; 

7.  the  redshift,  taken  from  the  LQAC  (blank  if  not  defined); 


7.  Summary  and  perspectives 

The  LQRF  is  build  to  provide  an  optical  representation  of  the 
ICRF,  retaining  the  ICRS  concept  of  being  defined  by  extra- 
galactic  objects.  Starting  from  the  LQAC  entries,  quasars  were 


A.  H.  Andrei  et  al.:  The  LQRF  -  an  optical  representation  of  the  ICRS 


403 


Fig.  19.  Sky  density  of  the  LQRF.  The  counts  are  in  bins  of  10°.  The 
void  regions  shown  by  white  circles  are  in  the  galactic  plane.  The  dens¬ 
est  patches  lie  on  the  SDSS  region. 


Fig.  22.  Distribution  on  right  ascension  of  the  A6  LQRF  minus 
ICRF  residuals.  There  is  no  systematics  above  4  mas  and  the  median 
is  105.3  mas. 


Fig.  20.  Vectorial  distribution  of  the  systematic  deviations  (north  up, 
east  right)  of  the  LQRF  to  the  ICRF  within  bins  of  10°.  The  average 
value  is  32.7  mas,  graphically  represented  by  the  off  map  arrow.  All 
large  values  however  befall  to  bins  with  5  or  6  sources  only. 
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Fig.  21.  Distribution  on  right  ascension  of  the  Atrcos^  LQRF  minus 
ICRF  residuals.  There  is  no  systematics  above  4  mas  and  the  median  is 
108.3  mas. 


Fig.  23.  Distribution  on  declination  of  the  Acr  cos  6  LQRF  minus  ICRF 
residuals.  There  is  no  systematics  above  4  mas  and  the  median  is 
108.3  mas. 


Fig.  24.  Distribution  on  declination  of  the  Ad  LQRF  minus  ICRF  resid¬ 
uals.  There  is  no  systematics  above  4  mas  and  the  median  is  105.3  mas. 


indentified  in  the  USNO  Bl.O,  GSC2.3,  and  SDSS  DR5  catalogs. 
The  positions  there  assigned  were  ensured  to  be  homogeneous 
with  respect  to  the  HCRF  by  applying  local  corrections  obtained 
in  small  neighborhods  around  the  quasars,  relative  to  the  UCAC 
and  2MASS  catalogs.  The  global  orientation  was  then  adjusted 


to  the  ICRF  and  the  zonal  departures  were  removed.  The  ICRF 
was  represented  by  the  ERF,  which  was  formed  by  accurate  long 
base  interferometry  radio  positions,  collected  from  the  ICRF- 
Ext2,  the  VCS6,  and  the  VLACalib,  and  selected  to  precision 
better  than  10  mas. 
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Fig.  25.  Example  extract  of  the  LQRF  catalog.  Its  columns  contain  the  elements  as  follows.  The  J2000  right  ascension  (h,  m,  s).  The  J2000 
declination  (d,  ").  The  external  right  ascension  error  (x  cos  S)  (mas).  The  external  declination  error  (mas).  The  number  of  radio  interferometry 

neighbor  positions  used  to  evaluate  the  external  errors.  The  R  magnitude,  taken  from  the  LQAC.  The  redshift,  taken  from  the  LQAC  (blank  if  not 
defined)  The  right  ascension  olTsets  to  the  radio  interferometry  position  (blank  if  there  is  no  radio  position)  (x  cos  S)  (mas).  The  declination  offsets 
to  the  radio  interferometry  position  (blank  if  there  is  no  radio  position)  (mas).  The  J2000  radio  interferometry  right  ascension  (h,  m,  s)  (blank  if 
there  is  no  radio  position).  The  J2000  radio  interferometry  declination  (h,  m,  s)  (blank  if  there  is  no  radio  position). 


The  final  LQRF  J2000  equatorial  cordinates  were  derived  by 
weighted  averages  of  the  input  catalog  positions,  once  locally, 
zonally  and  globally  corrected  as  indicated.  The  errors  assigned 
to  the  LQRF  coordinates  reflect  the  departure  from  the  ERF. 

The  Large  Quasar  Reference  Frame  (LQRF)  formed  in  this 
way  contains  100  165  objects,  of  which  2142  have  accurate  radio 
interferometric  positions.  The  comparison  between  the  LQRF 
positions  and  the  ERF  positions  indicates  that  the  overall  ori¬ 
entation  towards  the  ICRF,  as  represented  by  the  Euclidean  di¬ 
rection  cosines  is  Al  =  h-2.1  +  3.8  mas  A2  =  -0.9  +  3.5  mas 
A3  =  -2.6  +  3.4  mas,  with  zero  equatorial  bias  to  the  level  of 
2.9  mas.  The  average  offsets  to  the  ICRF  are  Aacosd  =  -1-2.7  + 
2.9  mas  and  Ad  =  -t-0.3  +  2.9  mas,  with  standard  deviations  of 
ctq,  =  134.3  mas  and  cr^  -  131.1  mas.  The  internal  errors  are 
described  well  by  a  Poisson  representation,  peaking  at  139  mas 
in  right  ascension  and  at  130  mas  in  declination.  The  LQRF  is 
planned  to  be  maintained  and  enlarged  in  the  future,  as  new  ver¬ 
sions  of  the  LQAC  and  other  quasar  input  catalogs  appear,  as 
well  as  to  incorporate  newer  versions  of  the  UCAC  stellar  refer¬ 
ence  frame. 
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